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ABSTRACT 


A  method  to  find  approximate  solution  of  an  idealized  model  of  the  blast 
wave  problem  is  developed,  the  gasdynamic  motion  is  initiated  by  an  instan¬ 
taneous  energy  release  from  a  point  source  (or  line,  or  plane  source) . 
Transformation  of  variables  is  Introduced  with  a  similarity  parameter  x  as  an 
independent  variable  along  with  a  quantity  less  sensitive  to  the  phenomenon. 

4 

This  reduces  the  fundamental  system  of  equations  of  the  problem  to  a  manage¬ 
able  form  for  the  approximation,  which  is  conveniently  performed  in  a  power 
series  expansion  in  the  variable  y  and  the  coefficients  of  the  expansion  can  be 
determined  successfully  from  the  systems  of  ordinary  differential  equations. 
The  method  has  been  applied  to  many  problems  of  the  similar  type,  with 
some  modifications  necessary  for  each  case,  to  blast  wave,  the  propagation 
of  the  blast  wave  in  the  non-uniform  medium,  exploding  wire  phenomenon, 
magnetohydrodynamics  cavitation  and  thunder. 


BLAST  WAVE  THEORY 


Alcira  Sakurai 
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1.  Introduction 

A  blast  wave  is  a  rather  common  phenomenon  usually  experienced  as  a 
"shock”  when  some  ex'^'osion  occurs,  and  the  phenomenon  itself  is  simply  a 
kind  of  disturbance  in  the  atmosphere  like  a  sound  wave.  The  characteristics 
of  a  blast  wave  are,  however,  quite  different  in  many  ways  from  those  of 
ordinars^  sound  waves.  Unlike  the  velocity  of  sound  c,  the  velocity  of  blast 
wave  U  is  not  constant  and  is  always  bigger  than  cj  usually  U  is  very 
large  near  the  source  of  the  explosion  and  decreases  very  quickly,  approaching 
the  sound  velocity  c  as  shown  in  Figure  1.  The  fact  Implies  also  the  energy 
dissipation  is  more  significant  in  a  blast  wave  than  in  a  sound  wave.  Secondly 
a  blast  wave  is  not  really  a  wave  of  periodic  type  like  ordinary'  sound,  but  con¬ 
sists  of  a  single  pulse  distinguished  by  the  presence  of  the  shock  wave  (The 
terminology  "shock  wave"  is  ustJd  here  to  indicate  the  front  surface,  not  the 
whole  pulse  region  which  we  refer  to  as  a  "blast  wave") .  At  the  front,  the 
pressure  p,  the  density  p,  and  so  forth,  jump  abruptly  from  their  values  at 
the  undisturbed  atmosphere  p^,  p^, ...  (  Figure  2) ,  Moreover,  its  wave  form 
changes  shape  in  the  course  of  propagation  in  quite  a  different  way  from  that  in 
the  sound  wave  case. 
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Figure  2 

Pressure  distribution  in  blast  wave 
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These  features  result  from  the  non-linearitv  of  the  phenomenon,  causing 
the  atmosphere  to  be  so  disturbed  that  the  pressure,  the  density,  and  so  forth, 
are  considerably  different  from  their  undisturbed  values  p^,  p^,  . . .  .  In  the 
sound  wave  case,  on  the  other  hand,  these  differences  are  small  and  the 
phenomenon  can  be  treated  much  more  simply  on  a  linearized  basis. 

Although  there  has  been  a  need  to  clarify  the  details  of  those  features  for 
many  years,  chiefly  with  the  obvious  aim  of  estimating  the  effects  of  explosives, 
the  attempts  at  theoretical  study  of  the  phenomenon  had  to  face  the  very  dif¬ 
ficult  problem  of  finding  the  solution  to  the  non-steady  flow  of  the  full  non¬ 
linear  hydrodynamic  equations,  satisfying  a  moving  boundary  condition  at  the 
shock  front. 

This  situation,  however,  has  been  eased  to  some  extent  since  the  Second 
World  War,  the  demands  of  which  obviously  stimulated  renewed  interest  in  the 
phenomenon.  As  a  consequence  a  large  amount  of  experimental  as  well  as 
theoretical  work  was  performed.  Considerable  progress  has  been  made  not  only 
in  this  particular  area,  but  also  in  related  fields.  The  investigation  has  since 
been  extended  to  quite  different  fields  by  applying  the  techniques  and  the  re¬ 
sults  proved  useful  in  the  original  problem. 

The  greatest  progress  was  made  probably  when  the  concept  of  similarity 
was  introduced,  which  simplified  the  problem  while  retaining  the  essential 
nature  of  non-linearity.  The  concept  itself  is  not  new,  but  has  been  familiar 
also  in  some  other  fields  of  fluid  dynamics,  such  as  boundary  layer  theory, 
conical  flow  theory  and  so  forth.  The  assumption  of  simileurity  causes  the  number 
of  independent  variables  to  be  decreased  and  thus  frequently  reduces  the 
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fundamental  partial  differential  equations  to  the  more  manageable  ordinary  dif¬ 
ferential  equations. 

The  blast  wave  solution  of  the  equations  of  gas  dynamics  was  found 
(Sedov,  1946,  Taylor,  1950)  in  the  form  of  similarity  solution.  However  the 
existence  of  a  group  of  similarity  solutions  of  progressive  wave  type  had  been 
known  before  in  connection  with  the  contracting  spherical  wave  problem. 

( Guderley,  1942) .  It  should  be  noted  that  the  similarity  in  the  blast  wave 
phenomena  does  not  hold  exactly,  but  is  valid  only  while  the  wave  is  strong 
enough  to  neglect  the  effect  of  the  ambient  atmospheric  pressure.  Since  it  is 
obvious  that  the  hydrodynamical  representation  of  the  phenomenon  does  not 
hold  very  near  the  explosion  source,  the  range  where  the  similarity  solution  is 
valid,  is  accordingly  limited  in  a  very  small  region  ( even  in  the  strong  blast 
wave  from  an  atomic  explosion,  it  is  reported  to  extend  between  20  and  180m 
(Taylor,  1950) ) ,  Sometimes  in  the  case  of  ordinary  explosions  there  is  no  such 
region.  Nevertheless,  the  concept  of  the  similarity  seems  to  be  very  useful 
since  the  phenomenon  preserves  the  similarity  in  some  extent  (  Sakurai,  1953, 
1954)  even  at  the  stage  where  the  front  shock  becomes  weaker. 

It  has  been  found  that  the  same  kind  of  similarity  solutions  could  be  obtained  for 
other  kinds  of  explosions  (Sedov, 1946,  Sakurai,  1953,1954,  Lin,  1954)  such  as  from  a  line 
source  and  a  plane  source,  while  the  ordinary  explosion  is  considered  as  that 
due  to  a  point  source.  Although  the  blast  wave  from  the  line  or  plane  source 
seems  to  be  rather  artificial  in  actual  explosions,  these  solutions  themselves 
proved  to  be  very  useful  in  many  other  different  fields,  especially  the  applica¬ 
tion  of  the  theory  to  the  hypersonic  flow  problem  (  Lin,  1954,  Lees  and  Kubota, 


5- 


1957)  where  the  general  principle  of  similitude  ( Hayes,  1947)  between  hyper¬ 
sonic  flow  and  blast  wave  behavior  had  been  knowm  to  hold.  Most  of  the  de¬ 
tails  on  these  earlier  developments  are  seen  in  such  books  as  Sedov  (1957), 
Hayes  and  Probstein  ( 1959) ,  and  the  more  extensive  treatment  on  the  subject 
given  in  the  recent  book  by  Korobeinikov  et  al.  ( 1961) . 

While  blast  wave  theory  itself  has  been  continually  improved,  many  im¬ 
portant  applications  of  the  theory  to  the  various  fields  of  research  have  also 
been  taking  place.  It  is  the  purpose  of  the  present  article  to  descirbe  some  of 
these  recent  developments.  The  fundamental  idea  of  blast  wave  theory  itself 
will  be  described  in  Chapter  2,  which  is  essentially  concerned  with  the  trans¬ 
formation  of  the  variables  -  both  dependent  and  independent  -  based  on  the 
concept  of  similarity,  and  reduces  the  fundamental  system  of  hydrodynamic 
equations  to  a  more  manageable  form,  although  they  are  still  non-linear  partial 
differential  equations. 

Ditterent  methods  of  finding  approximate  solutions  to  the  system  of  equa¬ 
tions  will  be  discussed  in  some  detail.  Various  applications  of  the  solutions 
thus  obtained  are  given  in  the  next  Chapter  3.  Firstly,  they  will  be  applied  to 
the  blast  v/ave  problem  itself  ( Section  3. 1) ;  this  proved  to  be  quite  successful, 
especially  in  describing  the  strong  blast  phenomenon  resulting  from  the  atomic 
explosion. 

Suppose  that  we  have  an  explosion  of  very  large  scale,  as  in  astrophysical 
phenomena.  It  would  then  be  necessary  to  modify  the  theory  so  as  to  include 
the  e^cts  of  non-uniformity  in  initial  frelds  caused,  for  example,  by  gravity, 

(  Section  3. 2).  [Applications  to  the  hypersonic  flow  problems  mentioned 
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previously  are  of  course  indispensable,  but  the  subject  may  be  omitted  here 
since  this  should  be  treated  more  extensively  in  an  independent  article.] 

The  electrical  disintegration  of  fine  wire  known  as  EWP  ( Exploding  Wire 
Phenomenon)  has  been  attracting  considerable  attention  in  connection  with 
various  fields  (see,Chace  and  Moore,  1959,  1962,  1964).  The  exploding  fine 
wire  may  be  regarded  as  an  example  of  the  line  source  and  may  produce  a  sim¬ 
ilar  situation  to  that  expected  from  the  blast  wave  solution  above  (  Section 
3.  3)  . 

In  the  vast  field  of  magnetohydrodynamics  (  Section  3.  4) ,  there  are  many 
problems  on  magnetohydrodynamic  shock  waves  which  show  blast  wave-like 
characteristics.  Since  extra  terms  are  needed  to  describe  the  effects  of  the 
magnetic  field,  the  genuine  blast  wave  theory  above  is  usually  not  applied 
directly  to  these  phenomena,  but  needs  to  be  modified.  Some  of  the  problems 
are,  however,  very  similar  to  those  of  ordinary  blast  wave  type  and  need  only 
a  slight  modification,  introducing  the  magnetic  pressure.  There  is  an  interest¬ 
ing  application  of  the  theory  to  the  problem  involving  singularity.  This  is 
given  in  Section  3. 5  in  connection  with  the  problem  of  the  collapse  of  an 
empty  cavity  in  water.  The  procedure  to  this  case  is  not  straightforward  be¬ 
cause  of  the  singularity  but  can  be  modified  by  introducing  Lighthill’  s  tech¬ 
nique  ( 1949) .  Quite  recently  application  of  the  blast  wave  theory  from  the  line 
source  has  appeared  in  the  problem  of  thunder.  Observations  of  this  have 
revealed  that  the  duration  between  the  lightning  and  the  succeeding  thunder,  its 
wave  form,  and  so  forth,  are  very  different  from  the  results  expected  from  the 
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theory  of  sound.  A  more  sophisticated  way  of  treating  the  phenomenon  is  re¬ 
quired  to  take  into  account  its  non-linearity,  and  to  regard  it  as  a  blast  wave 
from  the  lightning  as  a  line  source  (  Section  3.  6) . 

Apart  from  the  applications  mentioned,  blast  wave  theory  is  mathematically 
a  kind  of  transformation  of  variables  to  another  set  where  some  of  the  inde¬ 
pendent  variables  are  not  sensitive.  The  extension  of  the  theory  along  the 
line  as  a  technique  to  handle  the  complicated  problems  especially  of  non-linear 
nature,  could  be  very  useful  and  the  development  of  the  technique  itself  may  be 
an  interesting  problem  for  the  future. 

§2.  Blast  Wave  Theory 
2. 1  Fundamental  Equations 

Suppose  that  we  have  an  explosion,  following  which  there  may  exist  for  a 
while  a  very  small  region  filled  with  hot  matter  at  high  pressure,  which  starts 
to  expand  outwards  with  its  front  headed  by  shock  v/ave.  The  process  usually 
takes  place  in  a  very  short  time  after  which  an  advancing  shock  wave  develops, 
which  is  continuously  absorbing  the  ambient  air  into  the  blast  wave.  Although 
some  of  the  explosive  products  may  still  remain  near  the  center,  the  amount  of  the  air 
absorbed  increases  with  time,  and  the  later  behavior  ofthe  blast  wave  may  well  be 
represented  by  the  following  model  of  the  shock  wave  at  the  front  and  the  purely  gas- 
dynamic  treatment  of  the  motion  of  the  air  inside,  which  is  assumed  to  bear  ideal  and 
non -viscous  adiabatic  exponent  y  • 

The  three  types  of  blast  wave  (  spherical,  cylindrical  and  plane)  all  have 

features  in  common  and  may  be  conveniently  treated  in  the  following 
unified  manner.  The  equations  of  motion,  continuity  and  energy 
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of  the  gas  behind  a  blast  wave  of  one  of  these  types  may  be  written 


Du 

Dt 


1  iP 

p  ar  » 


(2. 1) 
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where  u  is  the  particle  velocityi  p  is  the  pressure,  p  is  the  density  and 
u,  p,  p  are  functions  of  the  Eulerian  coordinate  r  ( measured  from  the  center) 
and  the  time  t  ( measured  from  the  instant  of  explosion)  •  The  coefficient 
«t  has  the  values: 

a  =  0  for  a  plane  wave, 

«  =  1  for  a  cylindrical  wave, 

«  =  2  for  a  spherical  wave, 
and  the  expression  D/Dt  denotes 

Dt  “  at  ^  ar  • 


Equation  (  2.  3)  is  conveniei  ly  transformed  into 


Dt 


.du  ,  au  * 


(2.4) 


where  we  have  used  the  equation  (  2. 2) . 

The  position  of  the  shock  front  is  represented  by  R  ( measured  also  from 
the  center) ,  which  is  supposed  to  be  a  monotonically  increasing  function  of  t 
and  is  related  with  the  shock  velocity  U  by 


(2.  5) 
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At  the  shock  front  where  r  s  R,  the  quantities  such  as  u,  p,  p  measure 
suddenly  from  their  corresponding  values  in  the  atmosphere  ( presumed  at  rest) 
with  pressure  and  density  p^,  .  The  discontinuity  conditions  is  given  by 

the  Rankine-Hugoniot  relations  (  see  for  example,  Courant  and  Friedricks,  1948) 
which  for  the  present  purposes,  are  most  oonventiently  written 


'“>r=R=^ 


2  2 
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where  C  is  the  sound  velocity  given  by  C 
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(2.  7) 
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Now  we  have  three  equations  (2.1),  (  2.  2) ,  (2. 3)  and  the  boundary  conditions 
(  2. 6)  at  r  =  R  .  We  need  another  condition  to  determine  U  ( or  R)  as  a 
function  of  t  •  The  condition  may  be  postulated  in  various  ways  depending  on 
the  feature  of  the  explosion.  We  may  assume,  for  example,  a  small  region  of 
high  temperature  and  pressure  at  t  =  0,  from  which  blast  wave  is  started.  Some 
conditions  of  this  type  are  useful  especially,  for  the  purpose  of  finding  the 
solution  by  numerical  procedure,  where  all  the  complications  can  be  taken  into 
account  under  the  specified  circumstances. 

To  describe  the  general  features  of  the  phenomenon,  the  point  source  model 
( line  or  plane  source  for  or  =  0, 1)  has  been  considered  simple  and  appropriate 
for  the  purpose.  It  is  assumed  that  a  certain  amount  of  energy  is  released 
instantaneously  from  a  point  ( or  line,  or  plane) .  Since  the  total  amount  of 
energy  ( kinetic  and  thermal)  carried  by  the  blast  wave  is  constant  and  must 
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always  be  equal  to  the  energy  released,  the  assumption  of  the  instantaneous 
release  of  the  energy  is  conveniently  represented  in  the  following  way  ( Taylor, 


1950) . 


r  2  ’’■•’o.  a. 


or 


E  =  J  (r  pu  +  -*7)  r  dr - - 

«  '  2  v-r  v-1 


V“1  or+1 


(2.  8) 


where  £  is  related  to  the  released  energy  as 


Explosion  energy  per  unit  area  for  or  =  0  , 
( Explosion  energy  per  unit  line)  (  2ir)  ^  for 
( Explosion  energy)  (  4^)  ^  for  or  *  2  , 


a 


1  . 


(2.  9) 


Although  the  model  is  apparently  inadequate  to  represent  the  actual  situation  in 
the  very  early  stages,  where  all  kinds  of  complicatirns  are  involved,  it  should 
be  recalled  that  in  this  period,  the  gasdynamical  representation  of  the  phenom¬ 
enon  is  not  adequate  in  any  case.  The  advantage  of  using  the  model  is  two¬ 
fold.  Firstly,  the  model,  although  it  is  not  adequate  in  the  earlier  stage, 
becomes  more  and  more  accurate  at  the  later  stages  regardless  of  the  kind  of 
the  explosive  or  the  feature  of  the  explosion.  Secondly,  the  shock  front  ap¬ 
pears  from  the  beginning  in  this  model,  since  the  explosion  is  assumed  to  take 
place  instantaneously  at  a  point  ( or  line  or  plane) .  Thus  we  may  avoid  the 
nuisance  of  the  shock  formation  and  at  the  same  time  we  may  expect  a  kind  of 
similarity  in  the  flow  field.  Because  of  this  similarity,  the  whole  system  of 
equations  can  be  transformed  to  a  mathematically  feasible  system  and  thus  can 
be  solved  approximately,  which  will  be  seen  in  the  following  Section  2.  2. 
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It  should  be  noticed  that  since  an  instantaneous  source  is  assumed,  the 
second  shock  may  not  be  developed  in  this  model.  While  the  second  shock  or 
more  shocks  appear  in  actual  explosions  and  they  themselves  provide  interest¬ 
ing  effects,  their  features  depend  considerably  on  the  sequence  of  the  explosion 
and  are  well  understood  in  oonnecticm  with  specified  initial  conditions.  It 
seems  more  appropriate  to  use  direct  method  ( for  example.  Holt,  1956)  or 
numerical  method  ( for  example,  Goldstein  and  v.  Neumann,  1955,  Brode,  1959) 
for  studying  their  features  in  individual  cases. 

2.  2.  Blast  Wave  Transfcxtnation 

Let  us  introduce  new  independent  variables  ( x,  y)  defined  by 


2 


and  transform  the  dependent  variables  as 

u  =  Uf<x,y) ,  p  =  pQy“^  g(x,y) ,  p  =  pQh(x,y) ,  (  2.11) 

where  f,  g,  h  are  new  dependent  variables  to  be  found  <  Sakurai,  1953)  • 

The  introduction  of  x  is  most  important  and  is  based  on  the  similarity 
nature  of  the  one-dimensional  flow  given  by  Eqs.  (  2. 1) ,  (  2.  2) ,  and  (  2.  3) , 
for  which  a  group  of  similarity  solutions  depending  only  on  x  exists  (  See  for 
example,  Sedov,  1957) .  Although  none  of  these  solutions  satisfies  the  condi- 
tions(2. 6)  and  (2.8)  exactly,  there  is  a  .on  satisfying  them  approximately 
when  the  front  shock  is  strong  and  (  Pq/p)  j._j^  is  negligible.  Since  the  front 
shock  decays  very  quickly  even  in  the  stronger  explosions,  i*-  is  necessary  to 
take  into  account  the  effect  of  the  pressure  neglected  above,  for  which  we  need 
to  retain  independent  variable  other  than  x  •  The  variable  y  in  (  2. 10)  is 
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chosen  so  as  to  fit  the  conditions  (  2. 6)  and  (  2. 8) ,  which  are  now  written  as: 

«i»y)  •  . 


and 


^  9<i.y)  . 

» 


where  we  have  put 


®  Pq 


(  2. 12) 


(  2. 13) 


(  2. 14) 


The  transformations  (  2, 11)  are  thus  partly  chosen  to  fit  the  conditions  (  2. 6) 


and  (  2. 8) ,  and  at  the  same  time,  to  match  the  requirements  of  the  similar¬ 


ity  solution. 

Since  y  «C  ( P^/p)  for  a  strong  shock  (  See  Eq.  (  2. 6} ) ,  the  solution 
u  r—K 

tends  to  the  similarity  solution  as  y*0(R-*0)  .  Also  as  R  tends  to  <»,  y 
tends  to  1(  U-^C)  .  The  condition  (  2.12)  is  not  much  affected  by  the  value 
of  y  ,  which  varies  only  between  0  and  1  and  we  may  expect  a  similar  in¬ 
sensitivity  to  y  in  condition  (  2. 13)  also.  Thus  the  variable  y  is  expected 
to  have  little  effect  in  the  solution.  As  a  matter  of  fact,  we  can  make  other 
choices  of  the  variable  such  as  R  (  See  Section  3.  2  below) ,  which  are  expected 
to  be  insensitive  to  the  solution. 

By  introducing  (  2. 10)  and  (  2. 11) ,  the  fundamental  equations  (  2. 1) ,  (  2.  2) , 

(  2. 3)  are  now  transformed  into 


-13- 


< 


-Xg  +  ((-x) -^  +  Xy^.-vgC^  +  ^)  , 


(  2. 15) 


8x 


where 


X  .B  ^ 

^ “ y  dR  • 


( 2. 16) 


and  \  is  considered  as  a  fimction  of  y  only,  in  fact  R  is  represented  by 
Eq.  ( 2«13)  as  a  function  of  y  and  we  find,  by  differentiation  of  Eq.  (  2. 13) 


with  respect  to  y  , 


(2.17) 


vdiere  we  have  put 
1 


/  +-^)X*dX  «  J  . 


(2.18) 


Although  the  system  (  2. 15)  includes  the  integral  J  in  (  2. 18)  through  X. 
in  (  2, 17)  and  still  looks  formidable,  it  will  be  seen  in  the  next  section  that  the 
system  is  much  more  easier  to  treat  than  the  original  one  because  of  its  insensi¬ 
tiveness  to  y  . 

It  is  also  worthwhile  to  note  that  the  system  is  represented  by  three  non- 
dimensional  numerical  constants  «,  y  and  Rq,  in  which  R^  is  most  important 
and  is  related  to  the  scale  of  the  explosion. 

2.3.  Solution  in  a  Series  Expansion  in  y 
2.  3. 1  General  procedure 

It  is  readily  seen  that  if  we  put  y  s  o  in  Eq.  <  2. 15)  i'i  is  reduced  to  a 
system  of  ordinary  differeitfial  equations  where  \  is  constant  and  equal  to  «+l , 
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while  the  oondition  (  2. 12)  for  this  case  provides  the  boundary  values  at  x  s  1 
and  Eq.  ( 2. 13)  gives  the  relation  between  y  and  R  ( or  U  and  R)  as 
yOCR  .  These  are  all  the  immediate  results  of  the  similarity  solution  of  the 
intense  explosion.  Now  it  will  be  shown  for  general  y  that  if  we  express 
fi  g,  h,  in  power  series  of  y  as. 


< 


V. 


f  =  f<">+yf<'>.y2f<2)^...  , 

g  =  g^  +  yg^  g^  .  , 

h.h<">+yh<ny2h<2)^...  , 


(  2. 19) 


where  f^  ( i  =  0,  1,  2,  • . . )  are  all  assumed  as  functions  of 

X  only,  then  these  f^  ,  g^  ( i  s  1,  2,  •  •  • )  can  be  determined  suc¬ 
cessively  starting  fix)m  the  similarity  solution  given  by  f^  ,  h^ 

(  Sakurai,  1953,  1954) .  It  should  be  noticed  that  this  is  not  always  possible  in 
a  non-linear  system,  but  often  quantities  for  i  =  1  or  more  appear  in  the  first 
approximation,  which  makes  the  successive  approximation  procedure  impossible. 
Using  the  expression  ( 2. 19) ,  the  integral  J  defined  in  Eq.  (  2. 18)  becomes 
formally, 


2 


J  =  Igt  1  +  0-1  y  +  v^y  +...)  , 


(  2.  20) 


where  we  have  put 
r  1 


2  (0) 


<i) 


=  h'».^’)x“dx.  ,2.E1) 


.(2) 


'  0  ^ 

+  *1  /(h<°)  f<l)^  +  2h<^M^)  f<°Mx®dx  , 
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With  use  of  e3q>ression  (  2.  20) ,  Eq.  ( 2. 13)  becomes 
R  sH-I 

^  1  "  :A  efl)  ( Y-1)^  ^  2  ^  +•••]> 


(  2.  22) 


vi^icdi  shows  the  y-R  relation  more  clearly.  X  in  Eq.  (  2. 17)  is  also  expanded 


as 


Xs(a  +  l)(l  +  Xjy+X2y  +...)  , 


(  2.  23) 


where  we  have  put 


1 


'l‘ Jq(»H)(Y-1) 


< 


>.,  =  2  <r  ,  , 


(  2,  24) 


Putting  Eqs.  (  2. 19)  and  (  2. 23)  into  Eq.  (  2. 15)  and  equating  the  coefficients  of 
the  terms  in  corresponding  powers  of  y,  we  get  the  following  systems  of  dif¬ 
ferential  equations. 

From  the  coefficients  of  the  y*  terms, 

Y  ^ 

,  (2.25) 


Y(f<  ^  f^  -  1  =  (X  -  f^  ^ 


(0)* 


(0) 


From  the  coefficients  of  the  y  terms, 
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(  2.  26) 


< 


+  f< '''  M  X  -  f(  0) )  £(“> '  )h<«  X,  !«»  h<“'  , 

h«»  f<«  ■  -  ( X  -  f<°> )  h<« '  =  -  ( h'®)  ■  h'o* )  f<‘>  -  ( f«»  ■  h'»  . 

Y  g*®*  f<»  '  -  <  X  -  £«»  )  g<»  •  =  -  (g«»  '  g<“> )  f‘»  -  f<°> ) 

+  (»  +  l)\jg<°>  , 

^  •  •  •  • 


Similarly  £q.  (  2.12)  yield  the  following  conditions: 


rf(0),n--S-  --2i 

f  '^1' 

(  2.  27) 

'  Y+l  ’  ®  ’  V+l’ 

h“>(i) 

(Y-l) 

(  2.  28) 

=  0,  g'^*(l)  =  0,  h’^* 

(1)  ; 

(  2.  29) 

(Y-l) 

The  first  system  of  non-linear  differential  equations  (  2,  25)  with  the  con¬ 
ditions  (  2.  27)  gives  the  similarity  solution  (given  by  Taylor  (1950)  and  Sedov 
( 1946) ) ,  and  the  value  of  the  integral  in  Eq.  (  2.  21)  is  determined  from 
these  f^^^ ,  g^  ,  h^  .  The  second  system  of  Equations  (  2.  26) ,  after  in¬ 
serting  the  values  of  f^  ,  g^  ,  h^  obtained  above,  becomes  a  system  of 
linear  inhomogeneous  differential  equations  for  f^  ,  g^  .  The  system 

also  includes  the  unlcnown  paiameter  which  is  related  to  by  the  re¬ 
lation  (  2.  24) .  The  solutions  f^  \  ^  of  Eq«  (  2.  26)  subject  to  the 
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boundary  condition  (  2.  28)  thus  include  ( or  a  ^  •  Inserting  these  solu¬ 
tions  in  the  right  side  of  the  second  equation  of  Eqs.  (  2. 21) ,  we  have  an  equa¬ 
tion  to  determine  .  With  this  value  of  X^ ,  >  h^  are  determined 

finally,  and  we  have  the  second  approximation  to  the  prob  lem  to  the  order  of 
y,  in  the  following  form, 

u  s  U  {f^  +  y  f^  +  O(y^)  }  , 

p  =  (Pq/v)  +  y  g^^^  +  O(y^) }  ,  ( 2. 30) 

p  =  pQ  +  y  h^  +  0(  y^)  }  , 

The  proced  're  in  the  third  and  further  approximations  is  the  same  in 
principle  as  in  the  second  approximation  above,  and  f^  ,  g^  ,  h^  ,  x^  for 
all  i  are  to  be  found  successively. 

While  it  is  naturally  not  easy  to  see  the  validity  of  the  solution  in  ex¬ 
panded  form,  some  verifications  will  be  obtained  simply  by  comparison  with 

corresponding  experimental  data  as  well  as  with  some  purely  numerical  solu- 

—1  X(cH-l) 

tions.  It  is  noted  that  singularities  of  the  form  y  in  p  or  y  ^  '  in  R 

(  see  Eq.  (  2.  30) )  are  removed  bef<xre  the  expansion,  where  only  the  parts  pre¬ 
sumably  insensitive  to  y  are  concerned.  It  is  not  likely,  however,  that  the 
esqpansion  is  good  as  well  up  to  y  =:  1,  where  all  quantities  are  subject  to 
another  kind  of  singularity  ( see  Section  2. 4  below) . 

In  fact  it  is  expected  that  R— »  as  y-^1  and  accordingly  firom  Eqs.  (2.13) 
and  (2.18),  J  must  approach  the  value  («+l)  Ny"1)  ^  as  y“*“l,  and  it  is 
hardly  to  be  expected  that  the  expansions  of  J  in  Y  behaves  in  chis  way.  The 
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difficulty  will  be  overcome  to  some  extent  in  the  following  Section  2. 4.  Another 
technical  way  of  improving  the  expansion  procedure  is  to  expand  Pq/p  instead 
of  p/pQ  as  in  Eg.  (2.19),  which  led  to  the  expansion  of  [l  +  2y/(Y-l)]  ^  in 
2y/(  y-l)  (  see  Eg.  (  2. 12) ) ,  Since  2/(  y  -  D  is  usually  rather  a  large  number  for 
ordinary  values  of  y  near  for  y  =  1.  4,  the  expansion  in  2y/ (  y-1)  is  expected  to 
be  poor  and  the  alternative  way  of  expanding  Pq/p  is  supposed  to  give  better  results. 


2.3.  2  First  approximation 

Eguations  (  2.  25)  are  known  to  have  two  intermediate  integrals  and  the 
resulting  first  order  differential  eguation  is  integrated  exactly.  Multiplying 
the  second  eguation  of  Egs.  (  2. 25)  by  (y-1)  and  subtracting  the  third  one 
yields, 


JO) 


g 


(0) 


which  is  readily  integrated  to  give 

(0).  ,(0)xu(0)"''^^^ 

g  ( x  -  f  )  h 


x 


y+1 '  y+1 ' 


Y 


> 


(  2.31) 


where  we  have  used  the  condition  ( 2.  27)  to  determine  the  integration  constant. 


The  second  integral  is,  from  the  energy  relation. 


=  0  , 


which,  in  the  present  case,  is  reduced  to 

.2 

^  -  X  i 

y(  Y 


A.  ( a!.°!  +  i  £(0^(0)  j  «  (0) ,  ^  i  f(o)  )] 

dxl  *  v(Y-l)  2  '  Vl  2 


=  0 


and  integrated  to  give 

f<°>  =x(« 


(0)^  11  n<0)  (0)^  (0) 

.  I  JLa — wi - 

^yy-l  (0)^'  2  ^y-lj^(O)^ 


(2,32) 
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(0)  (0) 

where  we  have  again  used  the  condition  (  2.  27) .  Eliminating  g  and  h 
from  Eqs.  (  2.  31) ,  (  2.  32)  and  the  third  equation  of  Eq.  (  2.  25) ,  we  get 

,>  (0)'  o  (0) 

2-V  f(0)  ""..fO)'  ^(0).,’-'  > 

which  is  evidently  integrated  by  putting  ^  /x  =  F(  say)  with  the  result 
..  .(0)  «(«+3)(v-l)  .  (a+l)(v-l).  .  .(0)  . 

*  +  2v4-1  ^ 


e^4-2e4-5)  v^  +  (  -3tt^4-  2e-f  l)v+  4(  e^-l) 


(  2y+«-1)  {(«+!)  Y“(  ) 


log  I X  - 


Of  +  3 


+  log  c  , 


(  2.  33) 


where  c  is  the  integration  constant  determined  by  Eq.  (  2.  27) ,  given  as 


log  c  =  2  log 


Y+1  2y+«-1 


1  JCL  jL  f  a^Za+S)  v^+(  ~3o^+2flH-l)  v4-4(g^-l) 
Y+1  (  2 Y+«-l)  {(  0+1)  Y-(  «-l) } 


log  |l»+i.r(  «-iJbi  I 

*  ‘  (0+3)  (  Y“l)  '  * 


It  is  also  known  that  this  kind  of  exact  solution  exists  for  more  a  general  type 
of  similarity  solutions  such  as  the  initial  distribution  obeying  a  power  law  in 
distance,  time  dependent  energy  supply,  and  so  on  (  for  more  details  for  ex¬ 
ample,  see  Korobeinikov  et  al»  ^  1961),  It  is  noticed  that  the  solution  (2.33) 
is  not  always  valid  but  needs  to  be  modified  for  some  values  of  y  and  ct  . 

For  example  y  =  2  is  singular  and  the  solution  (2.33)  is  modified  by  applying 
a  limiting  process  as  y*^  2  namely, 


r<0) 


-xF,  F 


2iWl__  -«-3 
2F  - 1  -  27 


(  2.  34) 


a  s  2,  y  =  7  is  also  singular  and  the  solution  is  simply  given  by 
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(  2.  35) 


.(0)  _x  (0)  7, 

f  -4,  9  =4- 


The  case  y  ~  2  is  interesting  because  plasma  under  uniform  magnetic  field  is 
expected  to  behave  somewhat  like  a  gas  with  y  =  2  (  Spitzer,  1956) ,  while 
y  =  7  is  an  appropriate  value  for  water.  Graphs  of  f^  x) ,  g^  x) ,  h^  ( x) 
in  the  typical  case  of  y  =  1.  4  (  air)  are  shown  in  Figure  3.  Little  divergence 


Solution  curves  of  f^  ( x),  g^  x) ,  h^  ( x) 
for  o  =  0,  1,  2  taking  y  =  1.  4 
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( 0) 

is  found  for  various  values  of  \  .  f  at  x  -  0  is  always  zero,  showing 

u(0,t)=0  (2.36) 

as  required,  a  condition  sometimes  used  to  determine  the  flow  field  for 
further  approximations  (  see  below  in  Section  2.  3.  4) . 

After  obtaining  the  f^  ,  g^  ,  h^  values  above,  the  y  -  R  relation 
( or  U-R  relation)  given  in  Eq.  (  2. 22)  is  obtained  in  the  first  approximation  by 
evaluating  the  integral  in  Eq.  ( 2. 21) .  Unfortunately,  Jq  cannot  be  inte¬ 
grated  analytically,  but  is  only  found  numerically.  Since  the  exact  expression 
for  f^  given  by  Eqs.  ( 2.  31) ,  (  2.  32) ,  (  2.  33)  are  rather  in¬ 

convenient  to  compute  it  is  sometimes  more  efficient  to  get  their  values  by 
direct  numerical  integration  of  the  original  system  (  2.  25)  starting  from  the 
values  given  in  Eq.  (  2.  27)  at  x  =  1  .  The  most  recent  numerical  data  on 
f(  9)  ^  gt  ^  are  given  by  Jones  ( 1962)  for  some  values  of  y  .  values 

computed  from  these  data  are  exact  to  three  digits  when  compared  with  the 
values  obtained  by  various  methods.  Some  of  the  values  are  reproduced  in 
Table  1.  The  divergence  between  these  values  indicates  the  difficulty  in  find¬ 
ing  the  numerical  values  of  f^  ,  h^  partly  because  of  their  rapid 

change  near  x  s  1  . 


2.3.  3  Further  approximations 

It  is  convenient  to  introduce  new  variables  ( i  =  1,  2, . . . ) 


by  the  relations. 


1  s(x-f  )9  ,g  =g  ,h  =h  x  •  (2.37) 


Equations  (  2. 26)  to  the  second  approximation  are  thus  reduced  to 
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23 


(  2.38) 


-(X -!<«);«■  + 


.(  0) 


X  0) '  _£2L ,, J 1)  .„<«,  + a+i jifL V 

2  x-f<®>  2x-f‘°>^' 


(  2.  39) 


<x-f‘°*)(-Y»>‘*’'  +  +‘‘*')  =(*+l){(V-l)«>*‘’  +  4'***  -  .  (2.40) 


where  we  have  used  Eq.  (  2.  25)  • 

The  condition  (  2. 28)  at  x  *  I  becomes 


/‘’(l)  =-^.  +*‘’(1)  x“’(l)  .  (2.41) 

Another  condition  (  2.  21) ,  eliminating  o*  by  Eq.  (  2.  24) ,  is  transformed  into 

/{vf'‘’>(x.f<“>)h“"/«4a!i’  ^(«+^f«»\(°»x‘«}x“dx  (2.42) 

0  Y-1  ^ 

'S^o'^{<»+i)(y-:)  ' 

It  is  readily  seen  that  the  combination  of  Eqs.  (  2.39) ,  {  2. 40)  written  symboli¬ 
cally  as  (2v-l)x  Eq.  (2.39)  -  2x  Eq.  (2.40)  gives 

(x-f<°>)  +  (2y-1)x‘‘*]  =(<»+l)[v*^>-2i|i*‘*+ (2y-1)x'^’+2Xj]  , 


which  is  integrated  to  give 


-  2+*  ‘’  +  (2y-1)  x*  +  2Xj  ■=  (2Xj-  ^  ^) 


(  2.  43) 


exp 


4 

(/ 


W1 


1  X-f 


(0) 


dx) 


where  the  integration  constant  is  determined  by  Eq.  (  2. 41) . 
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+  ^01  ■**  (Y-lXif+l)  ,  (2.45) 


(1) 


where  Qp  0^,  Iqj,  Iq2 


n _ X_  ,( 0)  J  0)  .  _i_  <  0) 

<*2*273*  ^  *Y-1«  « 


*01  ♦  *02  \  •  2  2^!3  /  *“”  ^  W  > «  -  *V**<**  • 

Although  Eq.  ( 2. 44)  are  linear  in  ,  4<^  ,and  ,  it  is  apparently  not  possible 
to  find  the  solution  analytically,  except  for  the  very  special  cases  as  «  s  2  , 

Y  «  7i  to  which  the  first  approximation  becomes  very  simple  ( see  Eq.  (  2.35) ) , 
and  Pp  Pg  are  reduced  to  simple  forms  proportional  only  to  x  S  thus 
I  *1'^ 9re  found  in  the  form  of  a  certain  power  of  x  .  ( Morawetz,  1954, 
see  also  Korobeinikov  ^ ,  1961)  •  Generally,  Eq.  ( 2. 44)  may  be  integrated 
numerically.  The  equations,  however  contedn  an  undetermined  parameter  X.^  , 
and  are  not  suitable  for  numerical  integration  in  their  original  form  in  Eq.  (  2. 44) , 
but  may  be  integrated  in  the  following  way.  First  we  split  <p^  into  two 


parts  as 


9 


<1)  .  Jl) 


^^>4.  V  X  V 


(  2.  46) 


where  4*2^^  independent  of  X.^ 


Substituting  Eq.  (  2. 46)  horn  Eqs.  ( 2. 44)  and  ( 2. 45) ,  we  get  the  following 
systems  of  equations. 
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> 


(  2.  48) 


+‘2“'  =Pi4”  +  ‘’2+2‘’*% 


Xl{/(  Qi4‘’+  Oz+y  >  **  +  ^02  - 


(  2.  49) 


=  +02^1«)x-dx-I„j+^ 


€H>1) 


Eqs.  (  2.  47) ,  { 2. 48)  may  be  integrated  numerically  from  the  initial  values  at 
X  s  1,  which  are  obtained  from  Eqs.  (  2. 41)  and  (  2.  46)  as, 


"0  • 


2y 


Using  the  values  of  ^2^^ »  4^2^^  thus  obtained,  may  be  de¬ 
termined,  from  Eq.  (  2. 49) .  Once  X^  is  determined,  tp^  can  be  eval¬ 
uated  from  Eq.  (  2.  46)  and  Is  given  by  Eq.  (  2.43) .  The  numerical  pro¬ 

cedure  has  been  performed  for  y  -  1*  ^  In  each  case  of  or  =  0,  1,  2,  and  the 
values  of  X^  obtained  are  -2. 138,  -1. 989,  -1. 918,  respectively,  showing  these 

Table  n,  -X^  Values 
from  Korobeinikov  &  Chushkin 


a/y 

1.1 

1.2 

1.3 

1.4 

5/3 

2 

3 

0 

2.3257 

2.  2437 

2. 1862 

2. 1433 

2.  0683 

2.0143 

1.9407 

1 

2.  0866 

2.  0424 

2.0092 

1.9836 

1. 9374 

1. 9043 

1.8632 

2 

2.  0010 

1. 9666 

1.9396 

1.9182 

1.8785 

1.8496 

1.8141 

values  are  almost  2.  Recently  the  fact  has  been  verified  also  for  different 
values  of  y  by  extensive  numerical  works  by  Korobeinikov  &  Chushkin  ( 1963) . 
These  values  obtained  by  them  are  listed  in  Table  II. 
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Functions 


for 


Y  *  !•  4  are  shown  graphically  in  Figure  4. 


Figure  4 

Solution  curves  of  ^  for 

•  *  0,  ly  2  taking  y  4 

Exactly  the  same  procedure  as  above  is  applicable  to  the  further  stages  of  ap¬ 
proximations.  Swigart  ( 1960)  carried  out  the  third  approximation  in  the  case 
y  •U4p  •  ml  •  The  value  of  obtained  is 
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^2  “  2.7373  (for  Y  =  l.  4,  or  =  1)  , 

while  the  functions  ^  ^  ^  or©  shown  in  Figure  5. 


Figure  5 

Solution  curves  of  <p^ 

for  Of  =  1  taking  y  =  1.  4.  (  Swigart,  1960) 

2.3,4  Other  Methods 

It  is  seen  in  Figure  5  that  varies  considerably  and  its  value  is  rather 

large  compared  with  .  This  indicates,  as  suggested  in  the  pro¬ 

ceeding  section  2.  3. 1. ,  that  better  results  will  be  obtained  by  expanding  in 
Pq/ p  instead  of  p/ p^  .  in  the  expansion  of  p^/p,  the  quantity  corresponding 


( Z) 

to  X  zero  at  X  s  1  and  may  behave  better  in  the  whole  range. 

The  quantities  as  well  as  ,  themselves  vary  slowly, 

but  f J 9  •••  slowly  varying.  In  fact  become 

infinite  as  0(  x  '  ) ,  while  proper  choice  of  satisfying  £q.  ( 2. 49) 

cancels  the  singularity  in  •  Thus  ^  remains  finite  at  x  =  0  .  It 

should  be  true  to  for  all  i  since  u(  0,  f)  s  0  (Eq.  (  2.36) )  •  The  fact  may 

be  used  to  determine  as  follows 

Um  ( l{  **  +  V  4*' )  =  0  ,  (  2.  50) 

lc-0  *  ‘  ' 

or 

lim  ( +  X  )  =  finite  . 
x-0  ^  ^  ^ 

In  fact  the  value  of  X^  determined  in  this  way  coincides  with  the  value  above 
in  Table  H.  In  actual  procedure,  expanded  in  series 

of  X  near  x  «  0  and  their  ooetficients  are  evaluated  by  fitting  these  expansion 
solutions  to  the  numerical  one  bom  x  a  1  •  After  these  coefficients  are  de- 
termined,  X^  is  determined  to  cancel  the  singularity.  In  any  case,  it  is  not 
convenient  to  find  which  become  bigger  as  x-^0  because  of  the 

singularity. 

An  alternative  procedure  of  integration,  by  means  of  which  the  difficulty 
may  be  avoided  is  as  follows:  Introduce  the  notation. 
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Equations  (  2. 44)  and  (  2.  45)  are  then  conveniently  written  in  the  following 
vector  form, 


X»=AX  +  a  +  b  , 


/  m  X  X  dx  ^  0+1)  (y-1)  ”  ^01  ^  ^0  "  ^02^  * 


Now  we  introduce  new  functions  Y  as  , 


Y* 


V 


I  ,  =  <p(x),  4)  s  4;(x)  , 


(  2.  51) 
(  2.  52) 


and  multiply  Eq.  (  2.  51)  by  Y  *  from  the  left  and  integrate  nt  from  zero  to  one 


to  obtain 


.1  1  .1 
j  Y*  AX  dx+J  Y*(a+X..  b)dx, 

A  A  A  A 


which  yields  after  an  integration  by  parts, 

1  ,  1 
/(  Y*'  +  r*A)x  dx  =  [  r*x]n-  /  Y*(a  +  )dx  ,  (2.53) 

0  0  ^ 

Suppose  that  we  chose  Y  ,  so  far  undetermined,  to  satisfy  the  equation, 

Y*’  +  Y*A  =  tn*x“  ,  (2.54) 

Eq.  (  2.  53) ,  by  use  of  Eq.  (  2.  52) ,  becomes 


(  2.  55) 


Now  we  choose  the  value  of  Y  at  x  =  0  as 


which  is  in  detail 


(<P  5  )^_Q  =  0  • 
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since  \  \\i^  are  supposed  to  be  finite  at  x  =  0,  we  may  simply  put 

5(0)  *5(0)  *0  .  (2.56) 

to  satisfy  the  condition. 

Equation  ( 2.  54)  does  not  include  and  may  be  integrated  in  a  straightforward 
way  from  the  initial  values  (  2.  56)  at  x  *  0  .  The  solution  Y  thus  obtained  is 
used  to  estimate  the  integral  and  the  value  of  [T*X]„  ,  in  Eq.  (  2.  55) ,  from 
A  Lc*.  is  determined  as, 

S  '  +  'oi  -  <Io  -  l02  A*h«ix)-‘  . 

Once  is  evaluated,  Eq.  (  2. 51)  is  now  integrated  from  the  initial  values  in 

Eq.  (  2«  41)  at  x  =  1  .  The  process  seems  to  be  much  easier  than  that  for  Eqs* 

(  2*  47) ,  (  2.  48) ,  since  are  expected  to  vary  little  in  the  whole 

range  of  0<  x  <  1  .  The  process  is  essentially  the  same  to  find  X^,  X^,  ... 

(2)  (2) 

as  well  as  ^  ,  4<  ,  • . .  in  the  further  approximations. 

2. 4  Approximations  vedid  in  the  whole  range  of  y 

Although  the  solution  in  a  series  ir  y  may  be  good  near  y-^O,  it  is  certainly 
not  valid  near  y^l  .  It  is  also  noticed  that  the  expi;ession  for  the  y  -R  relation  given 
in  Eq.  ( 2. 22)  truncated  at  a  certain  power  of  y,  becomes  singular  at  a  particular  value 
of  y  .  In  the  second  approximation,  R  -  «  at  y  »  -l/X^,  where  R  is  supposed  to  be  still 
finite.  In  the  slightly  different  expansion  procedure  of  expanding  1/X  instead 
of  X  (  Korobeinikov  et  al. .  1961),  namely 

(• +1)A  =  1  -  X.jy  +  +  ...  , 

the  y  -  R  relation  are  obtained  from  Eqs.  (  2. 16) ,  (  2.  22)  as. 
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*  (  rJI^VJq)  y  exp  {-  X^y  +  “  +  . . .  }  , 

or  ( 2.  57) 

Rq  1  2 

y(-^)  =  Jq  exp  {X^y --X^  y  +  ...}  , 


where  the  singularity  at  y  =  -  l/x^  for  the  second  approximation  above  may  be 
avoided,  since  an  exponential  factor  remains  positive  all  the  time.  Neverthe¬ 
less,  it  can  not  be  valid  near  y^^l,  where  R  should  be  infinite,  while  R  in 
Eq.  (  2. 57)  terminated  after  a  finite  number  of  terms  remains  finite  at  y  =  1  • 
An  attempt  to  find  an  approximate  solution  valid  in  the  whole  range  of  y 
(  Sakurai,  1959)  will  be  given  in  the  following.  Let  us  first  make  two  assump¬ 


tions: 


8h 


i)  f«Cx  and  ii)  Xy  /dy  is  negligible  in  the  second  equation  of  Eq.  (2.15), 

(  f-x)  +  Xy  =  -h  (~  +~)  . 

The  assumption  ( i)  may  be  seen  in  Figure  9  below  in  Section  3,1.  3,  where 
u/U  =  f  for  the  second  approximation  for  or  =  2,  y  =  !•  4  is  shown  for  various 
values  of  y  and  ail  curves  are  actually  very  similar  to  straight  lines  through 
the  origin  x  s  0  .  The  second  assumption  is  based  on  the  following  facts:  the 
term  becomes  small  near  y  =  0,  since  y  enters  in  as  a  factor,  it  is  also  small 
at  y  1  where  X  becomes  small.  Although  it  is  not  clear  that  the  term  remains 
nail  in  the  intermediate  values  of  y,  and  both  assumptions  are  based  on  rather 
vague  reasons,  they  simplify  the  whole  situation  remarkably  and  thus  make  it 
>sslble  to  provide  an  approximate  formula  valid  in  the  whole  range  of  y  . 


From  the  first  assumption  we  get  simply 


(  2.  58) 
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where  is  the  function  of  y  only  and  given  by  the  condition  (  2. 12)  at  x  s  1  , 

2 


(  2.  59) 


By  use  of  Eq.  <  Z,  58) ,  the  second  equation  of  Eq.  (  2. 15) ,  neglecting  the  term 
8hy 

Xy  /By  according  to  the  assumption  (il)  becomes 

1  Wi 
h  dx 


1  8h  ^0  1 


-  f  X  * 
-  ^0  ^ 


which  is  readily  integrated  to  give 


.  .  m 

h  =  hgX  , 


( 2. 60) 


where  we  have  used  the  condition  ( 2. 12)  and  put 


Y-1  Y-1 


Y-1‘ 


m  = 


( «+l)  f( 

1  -  f. 


(  2.  61) 


With  use  of  f,  h  given  in  Eqs.  (  2. 58) ,  (  2. 60) ,  the  remaining  function  g  may 
be  determined  from  the  first  equation  of  Eq.  (  2. 15) ,  which  is  reduced  to 

a  yh^  ii^i  l+-\-  f^)  +  Xy}x  , 


and  integrated  to  give 


m+2  . 

g  a  A(x  -  1)  +  g^  , 

where  the  condition  (  2. 12)  is  again  used  and  A,  g^  are  given  by 


(  2.  62) 


A  a 


Yh, 


m+2  '*0 


0  y+1  Y+1 


(  2.63) 


Since  f,  g,  h  are  all  given  simply  as  power  functions  of  x,  the  integration  of 
J  given  in  Eq.  (  2. 18)  is  readily  performed  and  we  get 


J 


eh„f:+-^)  + 


Oq-A 


m  +  3+«  '2  0*0'v-l'  •  (Y-l)(irt-l) 


which  is  more  conveniently  written  in  the  following  form: 

JaP+XQ  , 


(  2. 64) 
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where  P,  Q  are  known  functions  given  as 


P  = 


0  0 


1  -  f 


yVo 


m+3+a 


^  2  ^0  ^  (  m+2)  ( Y*l)  ^  ^  (  Y"l)  (  m  +  2  ^  0^ 


Vh, 


0  = 


( 


1 


i  Y*^)  ( ni+2)  '  m+  or+  3 


(  2.  65) 


Equation  (  2.  64)  is  used  to  eliminate  J  from  Eq.  (  2«  18) , 

to  obtain  a  differential  equation  to  determine  \  •  Practically  it  is  easier  to 
eliminate  X,  leading  an  equation  for  J  , 

•^  =  ((‘»+i)I-;^llII-y^3’^  .  (2-66) 

It  will  be  seen  from  Eq.  (2.66)  that  J  tends  to  (a  +  1)  ( yi)  when  y 

approaches  1  and  accordingly  \~*0  there.  Since  m“*0;  g^,  hQ"*l  Qs 
y— 1  (c.f.  Eqs.  (2.59),  (  2.61),  (2.63)),  we  get,  from  Eq.  (2.65), 


from  which  we  get  J~*(  «+l)  \  Y"1)  ^  ^  >  thus  X  is  guaranteed  to 

become  zero  at  y  =  1  .  Equation  (  2.  66)  may  be  integrated  numerically  starting 
from  y  =  0,  where  X  =  a  +  1  and  J  “  Jq  Sfiven  by  Eqs.  (  2. 64) ,  (  2. 65) .  Equa* 
tion  (  2.66)  is,  however,  singular  at  y  =  0,  since 


y 


dy 


V 

X(  Y-1) 


f 


and  J(l--=r^)"*0  as  y— 0  .  To  proceed  with  the  numerical  integration,  it  is 
then  necessary  to  find  the  values  of  (dj/dy)y_Q  or  (dX/dyl^^^  which  are  related 
in  the  following  way: 
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(  2.  67) 


Using  this  relation  as  well  as  Eq.  (  2. 66) ,  it  is  found 

It  is  interesting  to  note  the  value  of  (a+1)  ^  (dX/dy)^^^  ,  which  corresponds 
to  X^  in  Eq.  (2.23).  For  y  =  1*4  this  gives  -2.32,  '1.82,  -1.61  for  asO,  1, 2 
respectively,  and  these  should  be  compared  with  their  exact  values  -2. 138, 

-1. 989,  -1. 918  given  in  Section  2. 3.  3. 

With  these  values  given  by  Eq.  (  2. 67)  as  well  as  (J  )y^Q  ®  ^H-1  at  y  =  0  , 
Eq.  ( 2. 66)  was  integrated  numerically  for  two  cases  of  a  s  1, 2  and  X,  ob¬ 
tained  finally  from  Eq.  (  2. 64)  is  plotted  in  Figure  6,  where  X-curves  given  in 


X 


Figure  6 

Decay  curves  (  X  vs  y)  from  various  results  for  a  =  1,  2  . 
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series  of  y  determined  in  the  proceeding  section  in  their  first  and  second  ap¬ 
proximations,  and  the  curve  obtained  purely  numerically  by  Goldstine  and 
von  Neumann  ( 1950)  for  a  =  2  are  shown  for  comparison.  Knowing  J  as  a 
function  of  y,  the  R-y  relation  (or  U-R  relation)  is  given  directly  by  Eq.  (2. 13), 

,  V  cf 

^'1  R  j  ‘  u2  ’ 

or  (  2.  68) 

R.  (u!r 
•^0  \c2 

It  is  easy  to  see  in  this  expression  that  R  becomes  infinite  as  y  approaches 
1  because  J  goes  to  ( a+1)  \  y^i)  ^  there  as  mentioned  above.  Its  asymptotic 
behavior  near  y  =  1  is,  however,  a  little  different  from  the  exact  one,  which  is 
known  as. 


for 

o 

II 

a 

1-y -C 

< 

r-3/4 

for 

a  s  1 

( log  R)^ 

for 

a  =  2 

(Whitham,  1950,  in  Sedov,  1957). 

It  is  noted  that  the  expression  (  2.  68)  itself  is  exact  as  far  as  J  is  exact,  but 
J  in  this  approximation  is  not  precise  enough  at  y'^1  and  1  -  y  behaves  as 
R  instead.  Nevertheless,  Eq.  (2.68)  with  this  value  of  J  is  expected 

as  a  whole  to  give  a  rather  good  approximation  to  the  whole  range  in  U-R  re¬ 
lation.  Its  defect  at  y'^1  may  be  improved  in  a  similar  way  to  that  given  in 
Korobeinikov  et.  (1961)  by  modifying  J  locally  at  y'*^l  to  fit  the  exact 

asymptotic  behavior  there. 

In  spite  of  their  comparatively  simple  forms,  the  approximate  solutions  of 
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f,  q,  h  given  by  Eqs.  (  2. 58) ,  (  2. 60)  and  (  2. 62)  seem  to  represent  the  feature 
fairly  well  (  Sakurai,  1959)  •  Because  of  its  simplicity,  the  g  function  in 
Eq.  ( 2. 62)  was  used  to  analyze  the  wave  form  of  thunder,  which  needed  trans¬ 
formation  too  difficult  to  perform  with  the  more  elaborate  formula  (Thome,  1962) 

(  See,  Section  3. 6) .  This  approximate  solution  has  also  been  used  by  Jeanmaire 
( 1963)  to  study  the  flow  field  in  T-shock-tube.  Another  approximation  theory 
is  considered  by  OsiLia  ( 1960) .  He  introduced  a  concept  of  "quasi-similarity" 
based  on  the  insensitiveness  of  the  functions  f,  g,  h  with  the  variable  y  and 
assumed 

SL  ^  _L^  ^  1  ^‘0  ■ 

8y  '  fg  dy  *»  8y  '  Sq  dy  dy  '  dy  "  > 

where  f^,  g^,  h^  are  given  in  Eqs.  ( 2.  59) ,  (  2. 63) ,  and  (  2. 61) .  Using  the 
expressions,  the  fundamental  system  of  equation  (  2. 15)  are  reduced  to  a  system 
of  ordinary  differential  equations  where  y  comes  in  as  a  parameter.  Using  the 
above  assumption  again,  dj/dy  is  reduced  to  an  integration  including  f,  g,  h 
where  y  enters  as  a  parameter.  The  procedure  to  find  the  solution  for  any  shock 
strength  y  is  as  follows:  start  with  a  guess  to  X  and  integrate  the  system  of 
differential  equations  numerically  with  given  value  of  y,  the  solution  is  used 
to  evaluate  J,  dj/dy  in  Eq.  (  2.17)  to  get  a  new  X  value,  with  which  we  re¬ 
peat  the  same  process  until  we  get  the  same  value  as  the  proceeding  X  value. 

Actual  computations  were  carried  out  for  a  =  1  ( cylindrical  wave)  at  respective 
-I/2 

values  of  y  s  u/C  =  1.1,  1.2,  1.4,  1.6,  2,  3,  «  .  The  last  two  cases  have 

been  found  to  coincide  with  the  solution  in  series  of  y  given  in  Section  2.  3. 

The  solutions  were  used  to  compare  with  his  experimental  data  on  the  exploding 
wire  ( c.  f. )  Section  3.  3) . 
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§3. 


lication  to  Various  Problems 


3. 1  Blast  Wave 

Although  application  of  the  foregoing  theory  to  actual  blast  waves  should  be 
the  main  subject  of  this  article,  it  is  not  easy  to  give  a  comprehensive  survey 
of  it  with  full  use  of  experimental  data.  An  excellent  description  of  the  phenom¬ 
enon  is  given  in  the  book  "The  Effects  of  Nuclear  Weapons"  (Glasstone,  ed. ) 
(1V62)  and  a  comparison  of  the  theory  with  some  data  is  well  reviewed  in  the 
book  by  Korobeinikov  et  ,  ( 1961) .  Only  a  brief  description  will  be  given  here, 
by  displaying  the  relations  which  might  be  useful  for  making  a  comparison  with 
experimental  data. 

3.1.1  Characteristic  lengths.  Scaling  laws 

As  long  as  we  assume  the  point  source  model,  represented  by  the  condition 

(  2.  8) ,  the  only  characteristic  length  that  appears  in  the  entire  formulation  is 

Rq  given  by  Eq.  (  2. 14) ,  and  the  features  corresponding  to  the  different  are 

l/(  a+1) 

simply  obtained  by  scaling.  Since  R^  is  proportional  to  Ej  and  to 

the  spherical  case  is  roughly  proportional  to  the  weight  W  of  the  explosive  ( of 


the  same  kind,  of  course) ,  the  scaling  factor  is  expected  to  be  proportional  to 
.  This  is  usually  known  as  the  Scaling  law  for  blast  waves  and  is  known 


to  hold  for  fairly  wide  range  of  W  (or  E^)  .  (  See  for  example,  "Effects  of 
Nuclear  Weapons"  ( 1962)  pp.  127-146). 


It  is  noticed,  however,  that  R^  is  also  related  to  in  the  form 
— 1/(  0+1) 

RqACPq  ,  and  there  exists  another  scaling  law  concerning  the  ambient 

pressure  p^  .  Since  the  rate  of  explosion  energy  converting  into  blast  wave 
energy  depends  significantly  on  the  ambient  pressure  p^  (the  rate  usually 
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decreases  as  the  pressure  is  reduced),  this  scaling  law  can  not  be  so  accurate 
as  the  last  one  above,  but  nevertheless  we  may  expect  the  scale  of  the  blast 
wave  from  the  same  source  to  be  magnified  as  the  ambient  pressure  is  reduced. 
This  fact  may  be  useful  in  producing  a  situation  equivalent  to  a  strong  explosion 
in  a  laboratory  experiment  with  a  relatively  small  explosion,  where  the  ambient 
pressure  is  reduced.  The  technique  was  used  successfully  in  the  study  of  a 
blast  wave  from  an  exploding  wire  (  See  Section  3. 4  below) . 

It  is  stressed  that  the  role  of  is  not  merely  the  scaling  factor,  but 
characterizes  the  entire  phenomenon  and  much  information  can  be  obtained  by 
simply  estimating  the  value  of  by  the  formula  (  2. 14) .  The  length  R^  itself 
indicates  a  distance  in  which  U/C  falls  to  roughly  about  1.  8  for  a  =  2  (  See 
Figure  7  below)  and  this  corresponds  to  about  2.  6  of  the  value  of  the  over¬ 
pressure  ratio  at  the  shock  front  (defined  by  ( P“ Po)y-.|^/Po»  shock 

condition  in  Eq.  (  2. 6) ) .  These  figures  are  useful  when  estimating  a  rough 
picture  of  the  range  in  which  the  individual  blast  wave  is  effective. 

Apart  from  the  ideal  model  of  a  point  explosion,  there  are  many  other  charac¬ 
teristic  lengths  in  actual  circumstances,  such  as  the  dimensions  of  the  explosive 
and  so  on,  and  some  of  them  are  in  fact  important  especially  for  understanding 
some  details  of  the  phenomenon,  while  the  length  R^  gives  an  over  all  feature. 
Another  characteristic  length  not  directly  related  to  the  explosion  but  important 
for  applying  the  point  explosion  theory  to  the  actual  situation  is  a  distance  be¬ 
yond  which  the  point  source  assumption  is  valid.  It  is  suggested  by  Taylor 
( 1950)  for  the  spherical  case  that  the  distance  should  be  much  bigger  than  the 
radius  R  (  say)  of  a  sphere  of  ambient  air  equal  in  mass  to  that  of  the  explosive. 
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The  criterion  should  be  extended  also  to  other  the  cases  of  cylindrical  and 
plane  waves. 

3. 1,  2  Decay  of  the  Blast  Wave 

The  decaying  characteristics  of  a  blast  wave  are  usually  given  by  the  re¬ 
lation  between  the  shock  velocity  U  and  the  distance  R  such  as  illustrated 
in  Figtire  1.  With  use  of  the  values  Xp  . . ,  given  in  Section  2.  3.  3  , 
the  relation  between  the  velocity  U  and  the  distance  R  (  U-R  relation)  is 
given  by  Eq.  (  2.  30) .  In  Figure  7,  U-R  relations  for  y  =  1*  4  are  shown  to  the 


Figure  7 

Velocity  -  distance  curves  from  Eq.  (  2.30) 


-41- 


second  approximation.  For  a  s  1,  y  ^  have  the  third  approximation 

utilizing  the  value  of  7373  given  by  Swigart  ( 1960)  (  See  Section  2.  3.  3) , 

vdiich  was  used  to  improve  the  formula  by  Lees  and  Kubota  ( 1957)  for  the  hyper¬ 
sonic  flow  past  a  blunt-nosed  axisymmetric  body. 

Sometimes  it  is  more  convenient  to  consider  the  relation  between  the  peak 
pressure  P  (the  pressure  at  the  shock  front  x  s  1)  and  the  distance  R  (  P-R 
relation) ,  which  can  be  derived  by  eliminating  U  (or  y)  from  Eqs.  (  2.  30)  and 
(  2. 6)  to  give 


The  formula  in  the  second  approximation  for  y  =  1*  ^  becomes 

+  2. 07  for  flr  =  2  , 

+  2.16  for  asl  ,  (3.1) 

+  2.  33  for  or  =  0  . 

Equations(3. 1)  have  been  compared  favorably  with  some  experiemental  data 


r 

1.961 

m 

^r; 

(  I.33I 

it; 

0.69 

i  ni 

( Sakurai,  1954)  • 

A  P-R  curve  for  a  s  2  valid  for  a  wider  range  established  both  by  numerical  as 

well  as  experimental  data  may  be  found  in  books  such  as  "The  Effects  of  Nuclear 

Weapons"  (1962) .  Peak  pressure  P  goes  down  very  quickly  in  the  beginning 

for  small  R  up  to  about  R  =  R^,  where  (  P/Pq)  -  1  is  2.  6  as  noticed  above. 

Beyond  that,  the  rate  of  decrease  slows  down  and  a  long  distance  is  needed  for 

the  wave  to  attain  accoustic  characteristics.  In  fact,  the  over-pressure  ratio 

goes  down  to  0(  10  S  when  R  *  lOR^  approximately  and  the  value  of  R/Rq 

•2 

reducing  the  over-pressure  ratio  to  0(10  )  may  well  be  more  than  50  .  It 
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is  noticed  that  the  over-pressure  ratio  for  the  ordinary  sound  wave  is  usually  of 

—6  ”2 
the  order  of  10  ,  in  comparison  with  which  the  above  figure  of  10  is  still 

very  high  and  enough  to  give  people  a  kind  of  feeling  of  "shock",  although  it 

may  actually  give  no  destractive  effects  anymore.  The  region  where  people  feel 

"shock"  is  thus  expected  to  spread  over  quite  a  range  and  the  region  is  also 

roughly  corresponding  to  the  domain  where  the  non-linearity  of  the  phenomenon 

is  dominant. 

Mother  important  quantity  to  express  the  decay  feature  is  X.  given  by  Eq. 
(2.16)  as: 

2RdU 

y  dR  ‘  “  U  dR  * 

It  is  stressed  (  Sakurai,  1955b)  that  \  is  rather  useful  for  both  experimental 
and  theoretical  studies  of  the  phenomenon,  because  the  characteristics  of  dif¬ 
ferent  kinds  of  blast  wave,  corresponding  to  different  values  of  a  can  be  com¬ 
pared  in  a  single  diagram  as  shown  in  Figure  6,  where  numerical  as  well  as  ap¬ 
proximate  results  for  a  =  2  are  compared  with  some  experimental  data. 

For  a  more  direct  way  of  plotting  experiment  data,  R  -  t  diagram  is  used, 
in  which  the  theoretical  curve  can  be  obtained  by  integrating  Eq.  (  2.  30) ,  The 


relation  in  the  second  approximation  becomes 

R/Rn 


K,K=/  {■. 

*^0  0 


-(  a+1)  -  .  ,-l/2  , 


K  for  O'  =  0  ,1  are  found  explicitly  as. 


(3.  2) 


(  W  v(1+v)  -  log("/l+v  +  Vv" I  },  v  =  for  a  =  0  , 


R  .  2 


(  W  1- JqV  R^^  “  for  a  =  1  . 
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K  in  the  first  approximation  is  reduced  to 


a+3'Rp’  * 

and  the  formula  for  o  =  2  shows  a  good  agreement  vith  the  data  for  an  atomic 
explosion  in  the  range  of  R  =  20~  i85m.  (Taylor,  19H0) , 

Oshima  (1960,  1962)  thoroughly  checked  the  reJatiori  (  3.2)  in  the  case 
or  s  1,  Y  s  1. 4  with  his  experimental  data  on  blast  waves  from  wire  explosion 
( c.  £.  Section  3.  3  below) .  Figure  8  ( Oshima,  1960)  shows  the  range  where  the 


Figure  8 

Experimental  data  for  shock  arrival  time 
compared  with  various  theoretical  curves  for  a  =  1, 
Y  =  l. 4  (Oshima,  I960) 
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formula  is  valid,  displaying  his  approximate  solution  mentioned  in  Section  2,  4 


3.1,3  Flow  Field  inside  Blast  Wave 

The  velocity,  the  pressure  and  the  density  Inside  the  blast  wave  are  repre¬ 
sented  by  f,  g,  h  given  in  Eq,  (  2,11)  and  their  changing  behavior  with  respect 
to  time  can  be  seen  in  their  graphs  against  x  for  various  fixed  values  of  y 
(or  U/c).  In  Figures  9,  10,  and  11,  these  graphs  for  a  =  2  are  shown  for 
various  values  of  U/c  .  Those  corresponding  to  other  cases  of  a  =  0,  1  are 
very  similar  in  nature  to  the  graphs  above. 


Figure  9 


f(  X,  y) 


for  various  values  of 


=  u/c) 
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Figure  10 


It  is  noted  that  the  usual  means  of  measurements  of  the  pressure,  in  actual 
observation,  is  performed  at  a  fixed  point  and  expressed  as  a  function  of  time, 
while  the  above  curve  of  g  is  concerned  with  a  fixed  time.  Sometimes  it  is 
necessary  to  convert  the  above  curves  to  the  one  corresponding  to  the  fixed 
point.  (An  example  of  the  conversion  is  given  in  Section  3.  7) . 

3,  2  Non-uniform  Field 

It  is  sometimes  necessary  to  take  into  account  the  non -uniformity  of  the 
medium,  through  which  the  blast  wave  is  propagating.  Blast  wave  propagation 
in  a  star  is  a  typical  example  of  this  case,  where  the  Initial  ambient  pressure 
and  density  p^,  pQ  ahead  of  the  shock  wave  can  not  be  assumed  constant,  but 
are  non-uniform  because  of  gravitational  effects.  Even  in  the  case  of  the  ter¬ 
restrial  scale,  for  huge  blast  waves  such  as  those  produced  by  nuclear  explo¬ 
sions  it  is  necessary  to  take  into  account  the  fall  in  pressure  and  density  in  the 
atmosphere  as  the  altitude  is  increased. 

It  is  possible  to  modify  the  blast  wave  theory  for  a  non-uniform  medium  by 
simply  assuming  that  the  initial  density  ( and  the  pressure  as  well)  may  be  ex¬ 
panded  in  a  power  series  of  R  as 

pQ  “  ®1  ^  *  * 

where  a^(  i  =  0,1,  2, ... )  are  constants  assumed  to  be  known.  A  series  expan¬ 
sion  solution  similar  to  that  developed  in  Section  2.  3  is  then  modified  to  in¬ 
corporate  the  non -uniformity  effects  represented  by  a^,  a^,  ...  with  additional 
terms  in  the  series. 

The  approach  will  be  illustrated  below  in  the  case  of  a  spherical  shock  wave 
in  a  star  (  Sakurai,  1956) .  Suppose  the  initial  density  distribution  in  star  is 
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given  by  the  following  Emden*  s  equation  ( slightly  modified  from  the  original 
one,  see  for  example,  Chandrasekhar,  1939) : 


<3.  3) 


where  we  put 


pQ  R  ^0 

-V 


p  ,  p  are  the  values  of  p^,  p^  at  the  center  of  the  star,  and  I*  represents 
a  length  connected  with  the  dimensions  of  the  star  given  by  L  s  (  ,  )  , 


2vp^G 


while  G  is  the  constant  of  gravitation.  We  find  from  £q.  (3. 3) , 

D  «  1  -  A^z^+  ...  . 


(3.4) 


It  is  also  necessary  to  modify  the  fundamental  system  of  equations  to  include 
the  terms  from  the  gravitational  effects.  Thus  Eq.  ( 2.1)  is  modified  to  give: 


supplemented  by 


Du  1  8d  Gm 
Dt  “  “  p  8r  "  2  * 


8m  .2 
«  4iir  p  , 


while  £qs.  ( 2.  2) ,  (2. 4)  are  unchanged  and 

Dt  P '8r  ^  r  '  * 


Dt  ^'8r  ^  r  '  • 


(3.  5) 


The  condition  ( 2.8)  for  an  instantaneous  explosion  of  energy  is  also  modified 
to  include  gravitational  energy. 
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(3.6) 


r  /I  2^  1  £ 

J  (r  u  +  — 

2  Y-1  P 


Gin  t  .  2 . 
-j7*)  p4nr  dr 


R  p  Gm_  , 


where  is  defined  by 


l/?. 

The  shock  conditions  (  2.6)  are  the  same,  c  is  given  by  (  VPq/Pq^  ^ 
not  constant  in  the  present  case,  but  a  function  of  R,  which  will  be  expressed 
in  series  of  powers  of  z  using  Eq.  (3.  4)  •  We  have  another  condition  at  the 
shock  hront, 


( m)  .  (  3.  7) 

Now  to  express  the  solution  in  series  form,  it  is  not  possible  to  expand  it 

“2 

in  series  of  y(or  U  )  as  given  in  Section  2.  3,  since  U  in  the  present  case 
is  no  longer  necessarily  a  monotonically  decreasing  function  of  R,  but  may 
even  increase  as  a  result  of  decreasing  density  and  pressure  in  the  equilibrium 
state,  and  would  then  cease  to  be  single-valued.  Accordingly  we  can  not  use 

Eq.  (  3,  3)  instead.  As 
is  not  indispensable,  but 
used  alternatively.  Thus 


applied  in  Eq.  (  2. 11) , 


y  as  an  independent  variable  and  we  use  z  given  in 
is  noted  in  Section  2.  2,  the  choice  of  the  variable  y 
any  other  variables  insensitive  to  the  solution  can  be 
we  Introduce  independent  variables  x,  z  defined  by 


r  R 


and  transform  u,  p,  p,  m  in  a  similar  way  as  to  that 
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«  f(  X,  z) ,  ^  «  P(  z)  ^  g( X,  z) ,  «  D(  z)  h(  x,  z)  , 


m 

m 


B  M(  z)  n(  X,  z)  , 


(3.  8) 


where  we  put 


P  HI 

J^  =  P(Z),  ^.M(2),  ■n^=|np  . 

c  c 


Equations  (  3.  5)  are  then  transformed  into, 

r 

Ki-xy 


f  + 

+ 

ZA^Ssi  M 

dx 

az 

U  dz  * 

Yh  dx 

* 

8n  . 

,  2  D 

2  ^ 

Ox  * 

3z  ““ 
M 

x  h  , 

i  ah 
h  dx 

.  z  ah 
h  az 

z  dD 
Ddz  " 

.af  _  2f 
”  ax  "  X 

» 

iiSL 
g  dx 

+  liSL 
g  az 

A  a. 

■^Ddz  * 

2z  dU 

U  dz  “  " 

^9x  X  ' 

X 


while  the  condition  (3.6)  gives, 


^  Dz^  J  -  2y  MDz^  K  -  I  «  1  , 


where  we  have  put 

J-/(|hf^  +  ;^)x^dx,  K  =  /hnx^dx 

T  -  r*/-L  O  a2  MD,  2  ^ 

I~J^<^-2yA  — )z  dz  , 

The  conditions  (  2. 6)  and  (3.  7)  at  the  shock  front  (r  =  R)  go  to 

«1.Z)  , 

n(l,2)  .1  . 


(3.  9) 


(  3. 10) 


(3.11) 
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We  have  thus  reduced  the  problem  to  finding  the  solution  of  the  system  of 
equations  (3.  9) ,  (3. 10)  subject  to  the  boundary  condition  (  3. 11) .  In  these 
equations,  such  quantities  as  D,  P,  M,  C,  I  are  assumed  to  be  )cnown  in 
series  of  z  .  The  expression  for  D  is  given  in  Eq.  (  3,  4) ,  from  which  ex¬ 
pressions  for  P,  M,  C,  I  are  derived  as 


2  2  4  4  4 

P  =  1-  yA  2  +“yA  Z  +...  , 


M 


< 


(3.12) 


C?  1  t  i\  1^2 -2,3, 

-^  =  i-(y-1)a  z  +'Jq(y“1)  a  Z  +  ...  , 

3 

.  I=:^{i--^YA^Z^  +  ||(4v-3)A^z‘‘+...)  . 


The  solution  of  the  system  will  be  found  in  a  similar  manner  to  that  used  in 


Eq.  (  2. 19) ,  in  power  series  of  z  .  We  assume  first, 

r  4zf<«  4... . 

h=h*°’  +  zh'^>  +z^h*^’  +... 

.  .-(0)  (1)  ,Z„(Z)  , 

*rZ2^  *rZn  'Teae^ 


(  3.13) 


V. 

where  ,  g^ ( i  =  0,  1,  . . , )  are  all  functions  of  x  to  be  de¬ 
termined. 

Inserting  the  expression  (  3.13)  in  J,  K  in  Eq.  (  3. 10) ,  they  are  expanded  as 

J  =  Igd  +<»jZ  +  +...),  K  =  Kq(1  +  PjZ  +  +  ...),  (3.14) 

where  J^,  a^*  •  •  • »  ^1*  ^2’  *  *  *  constants  given  by  integrals  simi 

larto  Jq>  o'p  ^2*  '**  in  Eq.  (2.21).  With  use  of  Eqs.  (3.10),  (3.12) 

and  (3.14),  U  is  expressed  as 
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.2  ^0 


(  3. 15) 


Substituting  expressions  (  3. 12) ,  <  3. 13)  and  (  3. 15)  in  Eq.  (  3. 9) ,  and 
comparing  the  same  powers  of  z  on  both  sides,  we  may  get  systems  of  equa¬ 
tions  for  f^  ,  g^  ,  h^  ,  n^  similar  to  that  given  in  Eqs.  (  2.  25) ,  (  2*  26) , 
In  the  same  way  we  get  from  the  condition  (  3. 11) , 


9“”(1)  h‘®Nl)  =1  , 

=  g*  **<«  »  {«  «o  , 


(3.16) 


•  •  e  a 


s  -  —  T  g^  (1)  *  -  h^^^  (1)  e-  T  i^\ 

Y+1^0»®  U)  *  ^^^^^2^0,  rf'*^l)  *0  , 


Now  the  first  system  of  equations  with  its  boundary  condition  given  in 

Eq*  (  3*  16)  is  the  same  as  the  first  approximation  given  in  Section  2.  3.  2  except 

( 0) 

for  the  additional  equation  for  n  ( x) , 

n<°>'  . 

The  second  system  for  ,  g^  ,  h^  turns  out  to  give  simply, 

f  eg  sh  an  aa^aO. 

Because  of  this,  the  following  systems  are  simplified  and  the  fourth  system  of 
equations  for  f  ,  g'  ,  h'  ,  n'  '  and  ««  is  reduced  to  exactly  the  same 


as  that  for  the  second  step  to  the  approximation  given  in  Section  2.  3.  3  except 

(  3) 

the  equation  for  n  .  At  the  same  time  the  third  system  of  equations  result- 

2 

ing  from  the  term  in  z  is  reduced  to  a  similar  but  slightly  different  system, 
which  can  be  solved,  and  the  numerioal  procedure  to  determine 
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x) ,  x) ,  ( x)  and  is  the  same  as  that  developed  in  Section 

2 

2.  3,  3,  thevalue  of  so  determined  for  y  =  4  being  0. 182  A  .  With  this 

value  of  as  well  as  the  values  of  X.^  for  y  ^  4  given  in  Section  2.3, 

the  velocity  of  front  shock  U  given  in  Eq.  (  3. 15)  becomes, 

=  0.  596  (I  -  0.82  -  1.14  Z^  ...) 


^  =1.30(-^)"^/^{l+0.4lA^(^)^  +  0.  57(-^)^  +  ...}  . 

C  0  0  0 

(2)  (2)  (2)  (2) 

It  is  readily  seen  that  functions  f  ,  g  ,  h  ,  n  and  the  constant 

2 

to  be  determined  in  the  third  step  are  all  proportional  to  A  ,  which  repre¬ 
sents  the  effects  of  non -uniformity  in  the  initial  distribution  given  in  Eq.  (3,4) . 
Thus  the  second  order  term  in  z  in  the  series  expansion  gives  the  non- uniformity 
effect,  while  the  first  and  the  third  terms  give  the  ordinary  attenuation  effect. 

In  the  further  approximations  beyond  the  fourth  step,  the  two  effects  are  not 
separated  but  enter  in  combination. 

Many  other  problems  may  be  treated  in  a  similar  way  to  that  illustrated 
above  as  long  as  the  non-uniformity  in  the  medium  remains  within  reasonable 
bounds,  so  that  it  can  be  well  represented  by  a  series  expansion.  The  expansion 
is  not  necessarily  to  be  symmetric,  (  such  as  above,  where  it  is  a  function  of 
R  only) ,  but  may  have  terms  including  variables  other  than  R  ,  The  solution 
may  be  assumed  in  an  expansion  form  fitting  the  series  expansion  to  the  initial 
distribution.  A  good  example  of  the  procedure  is  seen  in  the  problem  of  a  huge 


terestrial  explosion  such  as  that  produced  by  nuclear  bombs,  in  which  the  effects 
of  the  density  decrease  upwards  niust  be  taken  into  account.  The  initial  density 


pQ  may  be  assumed  in  this  case  as: 
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where  4  is  taken  as  a  vertical  coordinate,  with  4  =  0  at  the  center  of  ex¬ 
plosion  and  the  solution  is  found  in  series  expansion  form  as 


'^u/U  =  f*  ( x)  +  Z  f*  ( x,6)  +  . . .  , 


v/U  =  2  /^'(x,e)  +  ...  , 

(  p/Pq)  (  ^  ^  =  9*  x)  +  Z  9*  {  x,  6)  +  . . .  , 

p/Pq  =  h‘®*(x)  +  Z  h‘^’(x,e)  +  ...  , 


where  6  is  the  angle  between  the  direction  4  and  the  radius  r  (rcosG  =  4) 


and, 


2  =  U  s 


dt 


while  the  shape  of  the  shock  front  is  given  by: 


R  =  R  (1  +  Z  +  ...)  . 

More  details  of  the  procedure  are  found  in  Korobeinikov  et.  ,  ( 1961)  with  a 
slightly  different  initial  distribution  of  density 

pQ  =  a©  +  a^  4*^,  n;  constant  . 

There  is  another  way  of  approach  to  the  problem.  This  is  tc  find  the  exact 
solutions  of  the  fundamental  equation  by  seeking  a  similarity  solution,  which 
satisfies  the  equation  for  specified  initial  density,  determined  from  the  consis¬ 
tency  of  the  existence  of  these  exact  solutions  (  see  Korobeinikov  et.  al.  1961, 
also  Korobeinikov  and  Karlikov,  1960) .  Since  these  solutions  are  usually  valid 
at  constant  Mach  number  at  the  shock  front,  the  effects  of  the  changing  Mach 
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number  can  be  obtained  by  the  perturbation  method  similar  to  those  used  above. 

A  simple  example  of  these  specified  initial  distributions  is  given  by  «<  R 
for  certain  values  of  w  depending  on  y  ,  which  provides  a  similarity  solution 
for  a  strong  shock  wave  (  See  for  example,  Sedov,  1957)  .  Perturbation  from  the 
similarity  solution  by  series  expansion  in  z  ( or  y)  has  been  treated  in  a  simi¬ 
lar  manner  to  that  used  above  by  several  authors  (most  of  the  works  are  in 
Korobeinikov  ,  1961;  also  see  Rogers,  1956) ,  Among  them,  there  is  an 

interesting  case  to  which  the  solution  of  the  second  approximation  can  be  found 

( 3c^  t"  1)  t  v(  1  —  c^) 

exactly  where  w  =  y  +  I - ’’  ^  Korobeinikov  and  Ryazanov,  1959) .  This 

is  actually  the  generalized  version  of  the  special  case,  o-  =  2,  y  =  7  mentioned 
in  Section  2.  3.  3. 

Now  the  method  described  is  essentially  a  perturbation  procedure  starting 
from  a  similarity  solution  and  the  method  can  be  applied  to  many  other  different 
problems.  The  problem  of  the  propagation  of  shock  wave  produced  by  non-uniform 
motion  of  a  piston  can  be  treated  in  the  same  way  above  to  find  the  effects  of 
the  non-uniform  motion.  The  flow  caused  by  a  uniform  piston  motion  (not 
necessarily  a  piston  in  ordinary  sense,  but  possibly  a  uniformly  expanding 
cylinder  or  sphere)  can  be  expressed  by  a  similarity  solution  and  the  effects 
due  to  the  non-uniform  pitson  motion  are  incorporated  in  the  perturbation  terms 
similar  to  those  above  (  c.  f.  Kochina  and  Mel'nikova,  1958,  i960).  The  problem 
is  also  related  to  the  hypersonic  flow  around  a  blunt  nosed  slender  body .  More 
complicated  application  of  the  method  is  seen  in  the  problem  of  shock  propa¬ 
gation  due  to  non-uniform  piston  motion  in  a  conducting  fluid  with  a  uniform 
magnetic  field.  In  this  problem,  two  effects,  from  non-uniform  piston  motion 
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and  the  non-uniform  medium  appear  at  the  same  time,  since  the  interaction  of 
the  flow  with  the  magnetic  field  distorts  the  symmetry  and  thus  induces  an 
effect  similar  to  that  from  a  non-uniform  medium.  The  problem,  though  compli¬ 
cated,  can  be  treated  by  the  same  method  of  perturbation  as  that  described 
above  (Ness  eU  ^.,1963). 

In  dealing  with  the  propagation  of  shock  wave  in  stars,  we  always  encounter 
the  problem  of  determining  the  behavior  of  the  shock  wave  at  the  surface  of  the 
star,  where  the  density  goes  to  zero  and  the  solution  becomes  singular,  and  it 
is  necessary  to  find  a  solution  fitting  this  singularity  locally  (Gandel'man  and 
Frank-Kamenetskii,  1956,  Sakurai,  1960) , 

3,  3  Exploding  Wire  Phenomenon 

Shock  waves  produced  by  electric  wire  explosions  have  been  employed  by 
many  research  workers  as  a  useful  tool  in  the  diagnosis  of  the  exploding  wire 
phenomenon  or  used  directly  for  various  practical  purposes  ( c.  f.  Chace  and 
Moore,  "Exploding  Wires",  Vol.  1,  1959;  Vol.  2,  1962;  Vol,  3,  1964) .  Since 
the  exploding  wire  is  required  as  a  line  source  producing  a  cylindrical  blast 
wave  similar  in  nature  to  the  one  discussed  above,  attempts  have  been  made  to 
utilize  blast  wave  theory  for  the  purpose  of  studying  the  phenomenon  ( c.  f.  Sec¬ 
tion  "Shock  Waves"  in  "Exploding  Wires"  above),  which  shows  different 
features  depending  on  the  purpose  of  the  study.  Simple  application  of  the  theory 
seems  inadequate  (Bennett,  1959)  and  only  give  rise  to  confusion,  suggesting 
more  careful  consideration  of  the  situation.  It  will  be  shown  below  that  there 
is  a  proper  region  for  the  theory  to  apply  and  observation  in  this  range  shows 
good  agreement  with  theoretical  results  ( Oshima,  1960,  1962),  These  two 
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different  features  will  be  given  in  the  following. 


3,  3. 1  Applicability  of  Blast  Wave  Theory  to  Exploding  Wire  Phenomenon 

Since  the  shock  wave  produced  by  the  exploding  wire  was  noticed  to  be  a 
useful  tool  to  clarify  the  exploding  wire  phenomenon,  various  measurements 
have  been  performed  to  pick  up  the  shock  wave  by  different  experimental  tech¬ 
niques  ( c.  f.  "Explodi  Wires"  above).  It  has  been  found  that  measured  val- 

2 

ues  of  distance  and  arrival  time  of  the  shock  wave  plotted  in  an  R  -  t 

diagram  are  usually  almost  on  a  straight  line,  in  accordance  with  the  result  of 

the  theory  of  the  strong  cylindrical  blast  wave  {  c.  f.  Eq.  (3.2)),  and  the  amount 

2 

of  total  energy  of  the  wave  estimated  from  the  slope  of  the  line  in  R  -  t 
diagram  ( using  the  relation  (  2. 14) )  agrees  reasonably  will  with  the  correspond¬ 
ing  value  expected  from  the  other  estimates  ( Bennett,  1959b,  1962a,  Jones  and 
Gallet,  1962)  .  Nevertheless,  it  has  been  shown  by  Bennett  (1959),  that  there 

are  some  serious  deviations  between  the  more  detailed  results  of  the  experi- 

2 

mental  data  and  the  blast  wave  solution.  Slopes  of  the  line  in  R  -  t  diagram 
for  various  ambient  densities  under  the  same  energy  input  do  not  show  inverse 
square  root  dependency  on  the  density,  as  predicted  by  theory,  which  gives 

=  (  l/2)(  vEi/JoPg)  t 

(  Eq.  (3.  2) ) ,  (Oshima  ( 1960)  also  noticed  that  the  difference  in  ambient  pres¬ 
sure  hardly  affects  the  flow  field  inside  the  blast  wave  in  its  early  stages)  . 

Secondly,  the  careful  examination  of  the  data  reveals  that  vs  t  rather 

2 

than  R  vs  t  provides  a  better  fit. 

It  is  known  that  the  relation  R^*=<^t  for  n  different  from  2  holds  for  a 
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strong  shock  wave  from  a  source  with  energy  addition  varying  with  time  accord- 

i  -2 

ing  to  the  law  E^ct  which  is  usually  attained  by  an  expanding  piston 

motion,  (  c.  f.  for  example,  Korobeinikov  ,  1961) ,  Since  there  are  many 

reasons  to  believe  that  energy  is  being  supplied  in  the  course  of  growth  at 
disturbance  from  an  exploding  wire,  an  exponent  different  from  2  might  indicate 
a  similarity  flow  with  energy  addition.  However  it  is  not  clear  whether  the  ex¬ 
ploding  wire  behaves  as  an  expanding  piston  obeying  the  power  law  required  by 
the  similarity  solution  (  Bennett,  1959b) . 

Since  the  phenomenon  would  appear  to  be  a  far  more  complicated  one  than 
that  given  simply  by  the  line  source  theory,  it  seems  more  appropriate  to  find 
the  flow  field  numerically.  This  was  done  by  Rouse  { 1959)  based  on  a  model 
reproducing  the  anticipated  physical  situation  as  precisely  as  possible.  The 
numerical  results  were  compared  with  the  experimental  data  so  that  we  may 
justify  the  postulated  physical  picture.  Although  there  still  exists  a  limitation 
in  the  capacity  in  the  computing  machine,  it  may  be  possible  to  carry  out  such 
a  program  for  elaborate  models  of  the  type  considered  by  Bennett  ( 1962b,  1963) . 

The  above  remarks,  however,  do  not  necessarily  rule  out  the  possibility  of  the 

applicability  of  the  blast  wave  theory  to  exploding  wire  phenomenon,  provided  this 

is  used  in  the  proper  region.  It  is  recalled  that  the  blast  wave  theory  should  be  applied 

to  a  region  of  R  where  the  fundamental  assumption  of  instantaneous  line  ( or 

point)  source  is  satisfied.  The  criterion  is  roughly  given  by  R  »  R  (c.  f. 

Section  3. 1. 1)  where  R  is  the  radius  of  a  cylinder  of  air  equal  in  mass  to  that 

of  the  explosive.  Take  for  example  a  copper  wire  of  diameter,  0.1mm,  which 

gives  the  value  of  R  as  about  0.  4  cm;  accordingly  R  should  exceed  4  cm  or 
so  in  this  case. 
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Under  these  conditions  we  can  hardly  expect  to  find  any  relation  between 


the  experimental  data  for  R  <  4  cm  and  the  line  source  theory  to  this  case. 

2 

Now  there  is  another  limitation  due  to  the  relation  R  t  resulting  from  the 
strong-shock  assumption,  which  neglects  0{y)  terms  and  retains  only  the 
first  approximation  '  j.  f.  Section  2.  3) .  Suppose  we  impose  the  criterion  y  <0.1 
which  corresponds  roughly  to  R  <  0,  4  Rq  .  This  limitation  can  be  eased  by  use 
of  the  higher  approximations.  Even  with  these,  however,  it  is  usually  conven¬ 
ient  to  use  the  region  near  R  =  Rq  to  get  better  fit  because  of  the  nature  of  the 


R-t  curve.  The  magnitude  of  R^  in  the  cylindrical  case  is  given  as 

R  1/2 

^0  ■  2irp„’ 


where  E  is  the  explosion  energy  per  unit  length  of  the  line  source  and  is  ex¬ 
pected  to  be  5'^  100  joule/cm  in  the  typical  experiments  on  an  exploding  wire. 
The  values  of  Rq  for  this  range  of  energies  are  2.  8  12,  6  cm  at  Pq  =  1  atm. 

To  find  the  explosion  e  ^rgy  by  use  of  the  proportionality  relation  on  the 

2 

R-t  curve,  it  is  essential  to  use  the  experimental  data  in  the  region  for 

R  »  R  ,  R  ;<  0.  4  Rq  .  (  3. 17) 

It  often  happens,  however,  that  when  R  satisfies  the  condition  R  »  R  it  ex- 

2 

ceeds  0,  4  Rq  .  Therefore  the  region  where  the  proportionality  in  R  -  t  holds, 
may  not  exist  but  it  should  be  stressed  that  this  is  a  situation  where  the  higher 
approximations  should  be  used  to  estimate  the  energy  from  the  data  in  the  re¬ 
gion  of  R  >  0.  4  Rq  ,  Practically,  it  may  not  be  easy  to  find  data  correspond¬ 
ing  to  this  region,  because  of  limitations  in  the  capacity  of  ordinary  techniques 
to  pick  up  the  shock  wave,  since  these  usually  utilize  the  properties  prominent 
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in  a  stronger  shock  wave  such  as  its  luminosity.  There  is  also  a  more  funda¬ 
mental  difficulty  due  to  the  distortion  of  cylindrical  symmetry  in  the  flow  field 
at  the  distances  in  question.  Nevertheless,  the  measurements  satisfactory 
for  picking  up  the  weaker  shock  show  good  agreement  with  the  results  of  the 
blast  wave  theory  as  long  as  the  experiments  are  especially  designed  to  retain 
conditions  of  symmetry.  One  of  this  kind  of  experiments  will  be  given  in  the 
next  section. 

3,  4.  2  Cylindrical  Blast  Wave  from  Wire  Explosion 

Qshima  (I960,  1962)  performed  an  extensive  experiment  specially  directed 
forward  the  study  of  the  cylindrical  blast  wave  produced  by  a  wire  explosion . 
The  experiment  was  deliberately  designed  with  the  intention  of  studying  the  de¬ 
tails  of  the  blast  wave  itself.  Firstly,  wires  were  exploded  between  two  plates 
( actually  inside  the  chamber  of  a  thin  cylinder) ,  so  that  the  flow  field  could 
retain  its  cylindrical  symmetry  at  any  distance  from  the  wire.  Secondly,  the 
explosions  were  executed  under  reduced  pressures  (  See  Table  III  below)  which 
resulted  in  attaining  larger  values  for  and  made  it  possible  to  obtain  a 
region  in  R  satisfying  the  conditions  {  3. 17). 

Another  merit  of  having  larger  by  reducing  the  pressure  is 

to  have  more  accuracy  in  the  measurements  of  the  flow  field,  since  the  field 
is  in  fact  magnified  with  a  larger  characteristic  length.  It  is  noted  that  a  larger 
Rq  can  also  be  obtained  by  putting  in  more  electric  energy  into  discharge,  but 
increased  amounts  of  energy  would  more  likely  give  rise  to  other  troubles. 
Thirdly,  he  utilized  the  Mach-Zehnder  interferometer  to  measure  the  phenomenon. 
The  technique  is  familiar  in  experimental  studies  of  aerodynamics  and  can 
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provide  information  on  density  variation,  which  is  not  required  to  be  as  pro¬ 
nounced  in  other  methods  of  observation  such  as  these  making  use  of  lumin¬ 
osity,  and  the  measurements  could  be  carried  out  over  a  wider  range  in  R  , 
which  covers  not  only  the  range  where  the  series  expansion  solution  ( Section 
2,  3)  is  effective,  but  outwards  beyond  that  of  the  very  weak  shock  wave  re¬ 
gion.  Furthermore,  it  should  be  noticed  that  the  measurement  supplies  the 
whole  knowledge  on  the  density  distribution  in  the  entire  field  within  the  blast 
wave,  and  makes  it  possible  to  compare  the  distribution  with  the  one  obtained 
from  blast  wave  theory,  so  that  more  detailed  comparison  is  obtained  than  in 
the  distance-time  relation  (R-t  diagram) .  Bennett  { 1962)  also  used  the  inter¬ 
ferometer  but  he  utilized  it  to  study  the  exploding  wire  phenomenon  itself 
limited  to  the  range  of  very  small  R  . 


Various  cases  of  reduced  pressures  as  well  as  other  conditions  in  the  experiments 
carried  out  by  O  shima  are  shown  in  Table  III. 

The  energy  inputs  in  the  cases  of  I,  II,  III,  IV,  and  VIII,  are  almost  the 
same,  but  the  different  values  of  the  corresponding  reduced  pressures  are 
effective  in  producing  blast  waves  of  different  scales.  Density  distributions 
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measured  by  the  interferometer  in  these  different  cases  reveal  that  there  are 
some  distinct  differences  in  the  flow  features  at  the  respective  stages  of  the 
blast  wave.  Some  typical  data  for  density  distributions  at  the  different  stages 
are  reproduced  in  Figures  12,  13,  14,  from  Oshima' s  report  (1960) .  These 
should  be  compared  with  theoretical  h-curves  given  in  Figure  11  in  Section 
3.1.  For  strong  shock  such  as  in  Case  I  (where  the  Mach  number  is  larger  than 
about  2  but  not  so  large  as  to  violate  the  condition  R  «  R)the  flow  features 
as  shown  in  Figure  12  agree  closely  with  the  blast  wave  theory  almost  entirely 


Experimental  h  vs  x  distributions  for  case  I-l  ( Oshima,  1960) 
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Experimental  h  vs  x  distributions  for  case  IV-1  (Oshima,  1960) 
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in  the  region  of  x  ,  At  the  next  stage,  with  Mach  numbers  from  about  2  down 
to  about  1. 15  (  cases  II,  III)  (  See  Figure  13) ,  the  rarefaction  type  region  oc¬ 
cupies  a  considerable  part  of  the  center  of  the  blast  wave.  However,  the  flow 
field  of  the  remaining  part  agrees  well  with  the  point  source  theory  as  long  as 
we  use  solutions  valid  in  the  wider  range  of  R  (or  y)  given  in  Section  2.4) 
rather  than  those  given  in  the  strong-shock  theory.  (  Particularily  Oshima  used 
his  solution  given  in  Section  2.  4  for  the  purpose) .  At  the  later  stage  of  the 
weak  shock  region  such  as  in  cases  IV,  VIII  (the  Mach  number  is  lower  than 
1. 15  as  shown  in  Figure  14)  the  flow  pattern  is  very  much  complicated  by  the 
effects  of  the  succeeding  rarefaction  waves  and  the  secondary  shock  waves. 

Gk)od  agreements  between  the  theoretical  h-curve  and  the  observed  dens¬ 
ity  distribution  in  the  range  of  Mach  numbers  above,  suggest  that  we  may  de¬ 
termine  the  energy  value  by  utilizing  the  theoretical  relation  and  the  experi¬ 
mental  data  of  the  distance  R  and  the  arrival  time  t  .  Oshima  did  this  by 
plotting  the  data  in  a  diagram  of  R  -  R/Rq  in  the  following  way:  First, 

U  =  dR/dt  values  are  computed  from  the  R  -  t  data  so  that  we  get  an  experi¬ 
mental  relation  between  R  and  U/c  .  Secondly,  theoretical  R/Rq  values  can 
be  found  for  each  U/c  values  from  the  theoretical  U-R  relation  as  given  in 
Eq.  (2.22)  or  (Eq.  (2-57),  Eq.  (2-68)).  Suppose  we  plot  tiiCse  experimental 
R  values  and  theoretical  R/Rq  values  corresponding  to  the  same  U/c  values. 
We  may  expect  a  linearty  in  the  (R)  exp.  -  (R/R^ltheo.  diagram  as  far  as  the 
phenomenon  carries  the  characteristics  predicted  by  the  point  (line  in  this  case) 
explosion  theory.  The  diagram  given  by  Oshima,  using  the  second  approxima¬ 
tion  formula  Eq,  ( 2-  30)  as  the  theoretical  U-R  relation,  is  reproduced  in 
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Figure  15.  In  this  diagram  we  can  see  the  linearity  especially  for  the  cases  of 
stronger  shocks.  It  is  noted  that  the  deviation  from  the  linearity  in  the  figure 
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Relations  between  r/R^  and  R  based  on  the  second  approximation  in  Eq.  (2.30)  (Oshima,  1960) 


does  not  necessarily  mean  the  failure  of  the  point  source  theory  but  instead  it 
reveals  the  inadequacy  of  the  second  approximation  to  this  region.  The  slopes 
of  these  lines  give  the  values  of  Rq  for  respective  cases.  From  the  values  of 
Rq,  the  explosion  energy  ( or,  more  precisely,  the  energy  equivalent  to  produce 

the  blast  wave  generated  from  a  line  source)  can  be  obtained  from  the  formula 

2  — 

E  =  2ir  pQ  Rq  .  To  check  the  E  values  thus  obtained,  Oshima  estimated  the 

E  values  by  use  of  a  more  direct  method  of  evaluating  the  energy  integral  given 

in  Eq.  (  2-8) ,  This  method  used  the  observed  density  distribution  p  and  the 

values  of  u,  p,  estimated  from  the  p  values.  These  E  values  agree  very 

well  with  those  obtained  above. 

4  M  aqnetohvdrodvnamics 

Essential  features  of  the  blast  wave  phenomena  in  magnetohydrodynamics 
are  expected  to  be  similar  in  nature  to  those  discussed  above  for  non-conduct¬ 
ing  fluids  as  both  resulted  from  the  same  sort  of  non-linearity.  However  some 
modifications  are  needed  to  represent  the  interaction  between  velocity  and  mag¬ 
netic  fields.  In  general  two  different  approaches  are  used  depending  on  whether 
or  not  electrical  conductivitv  is  involved.  In  the  case  of  low  conductivity, 
which  is  usually  the  situation  in  laboratory  experiments,  the  modification  of 
the  theory  is  more  or  less  straightforward  and  the  perturbation  method  may  ef¬ 
fectively  be  used.  It  should  be  noted,  however,  that  the  procedure  is  not  uni- 
formally  valid,  but  is  valid  only  in  the  finite  region  of  the  space.  This  procedure 
will  be  given  in  the  subsection  3. 4. 1,  where  the  discussion  will  be  confined  to 
the  specified  type  of  flow  preserving  axial  symmetry.  This  is  rather  typical  in 
this  field  of  magnetohydrodynamics,  since  this  is  conveniently  realized  with 
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axial  or  azimuthal  magnetic  fields  in  connecting  with  recent  experiments  on 
pinch  effects,  exploding  wire,  etc. 

The  case  of  high  conductivity  may  well  be  approximated  by  the  limiting 
case  of  infinite  electrical  conductivity,  more  exactly,  infinite  magnetic 
Reynolds  number.  Since  the  magnetic  field  is  "frozen”  in  the  fluid  in  this  case, 
the  assumption  sim.plifies  the  situation.  However,  this  sometimes  brings  about 
confusion  leading  to  erroneous  conclusions  and  hence  we  must  proceed  with 
caution  along  this  line.  It  is  possible  to  have  a  completely  different  situation 
depending  on  whether  or  not  the  initial  magnetic  field  is  permeated  in  the  fluid . 
These  will  be  the  subject  of  the  Subsection  3.  4.  2. 


3. 4. 1  Blast  Wave  Phenomena  in  the  Fluids  of  Low  Conductivity 

We  consider  here  the  cylindrical  type  flow.  The  magnetic  and  electric 
fields  B,  E  consistent  with  and  preserving  the  symmetry  are  given  by 


B  =  (0,  Bq,  B^),  E  =(0,  Eq,  E^) 

where  B.,  E*;  B  ,  E  are  the  azimuthal  and  the  axial  components  and  are  as* 
o  t?  z'  z 

sumed  functions  of  r  and  t  only.  The  equations  of  motion,  continuity  and 


energy  given  in  Eqs.  (  2. 1) ,  (  2.  2) ,  (2.  4)  are  now  modified  to  give 


B  BB  B 

Du  ^  Iz  ox 

Dt  ”  ar  "  8r  ”■  pr  ar  ' ^  ■  > 


(3. 18) 


Dt  “ 


(  3. 19) 


Dp  /  \  .  V~1  1  a  ->  •> 


(  3.  20) 


where  a  is  the  electrical  conductivity,  is  the  magnetic  permeability  and 


they  are  assumed  constant.  These  are  supplemented  by  the  Maxwell  equations, 
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dB 

BE  e 

ar  ”  at  * 

(3.  21) 

SB 

7  a7^^^e^  '"“‘aT  ’ 

(3.  22) 

(  3.  23) 

3B 

IT  . 

(3.  24) 

Since  we  assume  <r  is  small,  the  ordinary  shock  condition  (  2.  6)  may  be  used 
at  the  front  r  =  R,  and  the  magnetic  and  electric  fields  By  E  are  continuous 
there,  thus 


(t),.R  =  to  (3.25) 

where  Bq,  Eq  are  the  values  of  B,  E  in  front  of  the  shock  wave  presumably 
known.  They  may  be  functions  of  R  . 

Suppose  the  magnitudes  of  the  given  field  B^,  Eq  are  small  or  of  the 
ordinary  order  in  a  sense  such  that  its  magnetic  pressure  is  smaller  than  or  at 
least  comparable  with  the  pressure  Pq,  then  it  may  then  be  expected  that  the 
interaction  effects  of  these  fields  to  the  flow  are  small.  In  fact  Eqs.  (  3.  23) , 
(  3.  24)  show  that  the  qualities 


7  dr 


ar 


are  small  when  o-  is  small  and  deviations  in  B  ,  from  their  unperturbed 

values  BqqjB^q  are  thus  expected  to  be  small.  The  first  approximation  to  these 

deviations  may  be  found  by  replacing  those  qualities  B,  £,  u,  (involved  in 

terms  multiplied  by  <r)  with  their  unperturbed  values  Bq,  Eq  srid  u  ( given  by 

the  solution  of  corresponding  non-conducting  motion) ,  Thus  we  get  from  Eqs. 
(3.18),  (3.20),  (  3.  23),  (  3.24), 
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^pp-V=<r(Y-l)  p'''[(E,o+«Bg„j"+(uB^O-E0g>"]  , 


r  8r  *^zO  ■*■  “  ®eo'  ’ 


ae 

^  dr  ■  ^  ®zO  "  ^eo^  * 


(3.  26) 


The  last  two  equations  of  Eqs.  ^3,  26)  give  the  disturbance  in  the  magnetic 

field  due  to  the  interaction  and  thus  we  have 
^  r 

r  Bg  =  Po-  /  (E,n  +  u  B^g)  rdr  +  R  B^^  , 


R 


zO 


eo^ 


80 


(  3.  27) 


®z  ■  ®z0  “  ^00^  ■*■  ®z0  * 

R 


The  rest  of  the  equations  in  Eq.  (3,  26) ,  supplemented  by  Eq,  (  3, 19) ,  provide 
a  system  to  determine  u,  p,  p  •  The  system  may  be  solved  in  a  fashion  simi¬ 
lar  to  Section  2,  3  by  expanding  the  solution  in  powers  of  y  .  After  which  the 
deviation  due  to  the  magnetic  field  comes  to  the  order  of  cr  in  their  respective 
terms  as  f^^^ ,  ,  h^  ( i  =  0,  1,  , . . )  . 

This  method  has  been  used  to  estimate  the  effects  of  the  applied  axial  mag¬ 
netic  field  to  the  cylindrical  blast  wave  produced  by  exploding  a  fine  wire 
{  Sakurai,  i962a,  Sakurai  and  Takao,  1963,  1964),  in  which  Bq,  Eq  are  assumed 
as 

®60  =  fo  =  ° 


and  B  given  in  Lq.  (  3.  27)  is  simplified  to, 
z  ' 

r 

B  =  Bf.  ( 1  +  Ij-cr  /  u  dr)  . 
z  0 


(3.  28) 
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Suppose  we  assume  the  unperturbed  flow  u  is  given  blast  wave  theory  in  Eq, 
(  2. 11)  •  Then  Eq.  (  3.  28^  has  the  more  explicit  form 

X 


B  = 
z 


( 1  +  jio-  UR  /f(  X,  y)  dx)  ,  x  <  1 
1 


(3.29) 


where  the  expressions  (  2. 10)  and  (  2. 11)  have  been  used.  Utilizing  the  series 
expansion  solution  ( Section  2.  3)  we  have 

f  ^f(O)  f(l)  ^ 


^\2 

UR 


=  In  +  •  ••]  » 


ahd  Eq.  (  3.  2)  gives 


s  1  +  R  ( /  f*  °'dx  +  y(  /f*  X.,  /  f*  dx)  +  . . .  ]  ,  X  <  1  , 

Bp  m  j  j  ‘  1 

where  we  have  defined  the  magnetic  Reynolds  number  R  as 

m 

o-^icR^ 


R  =  -7=-  . 

m  slj^ 

The  values  of  B^/B^  computed  to  the  first  term  (neglecting  0(y)  terms)  are 
plotted  In  Figure  16.  Electric  field  E^,  derived  from  Eq.  (3.23)  using  the  ex¬ 
pression  for  B  above  becomes  singular  at  x  =  0  as  0(  l/x)  .  This  shows  the 
z 

procedure  is  not  uniformly  valid,  as  noticed  above.  Perturbation  equations  for 
U)  P)  given  in  Eq.  (  3.  26)  become  in  this  case, 

Du  ^  -2 

PdT' ' 


Dp  ^  u»  .  .  _2  2. 

dT  '"WtaT  +<r(v-BpU  )  , 


(3.  30) 
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Figure  16 


Decrease  in  axial  magnetic  field  due  to 
its  interaction  with  cylindrical  blast  wave 


where  u  ( in  the  terms  multiplied  by  <r  )  is  supposed  to  be  given  by 
u  =  U  f(  X,  y)  =  U(  +  y  +  . . .  /  . 

It  is  convenient  to  introduce  new  dependent  functions  f,  g,  h  defined  by 
u=:Uf(x,y),  P  =  Prt  y  ^  g(x,y) ,  p  =  p  h(x,y)  , 


where  upon  Equations  (  3.  30)  are  transformed  into, 


-Xg  +  (f-x)-g+Xy|^  =  -Yi(||  +  i) 


0  2 

+  (v-l)cr^iUR— y  r  . 


By  expanding  f,  g,  h,  in  a  power  series  in  y  , 

f-f  +yf  +•••,  g-g  +y9  +«..,h 


{  3.  31) 


+yh<‘>  . 
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and  utilizing  the  relation 


o-jiUR  =  [1--  X^y  +  ...]  , 

we  get  systems  of  equations  from  which  to  determine  f^^^ ,  ,  h^  ,  i  =  0, 1,  2, . . . 

It  is  observed  that  the  terms  multiplied  by  <r  in  Eq.  (3,31)  all  originated  from 
the  first  power  of  y  in  their  power  series  expansion.  Therefore  the  deviations 
in  f ,  g,  h  from  the  unperturbed  f,  g,  h  do  not  appear  in  the  first  approximation 
which  neglects  terms  of  the  order  of  y  (  strong  shock  approximation) .  The  same 


kind  of  perturbation  process  may  be  applied  to  other  cases  such  as 

®eo=f>  b,„  =  e„  =  o  , 


zO 


(  3.  32) 


which  represents  a  situation  associated  with  the  problem  of  a  blast  wave  in¬ 
fluenced  by  a  magnetic  field  produced  by  an  electric  current  flowing  through  the 
center  axis.  It  is  not  necessary  that  the  current  be  constant  with  time  only  that 
it  is  a  slowly  varying  function.  More  complicated  situation  like 

®z0  "  ^0  "  ®  » 

may  also  be  treated  in  the  same  manner.  Note  that  this  case  is  important  to  the 
study  on  the  effects  of  the  applied  axial  magnetic  field  to  the  exploding  wire 
phenomena,  and  differs  firom  the  previous  case  as  is  not  zero.  This  com¬ 
ponent  should  be  considered  in  the  study  of  the  early  stage  of  the  phenomenon, 
since  the  electric  current  of  the  exploding  wire  produces  azimuthal  magnetic 
field  for  a  while. 

It  is  noted  that  the  case  given  by  Eq.  (  3.  32)  can  be  represented  by  a  similar¬ 
ity  solution  under  the  assumption  of  strong  shock  and  s  =  constant  ( Greenspan, 
1962,  see  also  Korobeinikov  and  Ryazanov,  1962) .  The  solution  has  the  following 
form, 
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(3.33) 


u  =  Uf(x),  p  =  pp^g(x),  p  =Poh(x), 


While  the  expressions  for  u,  p,  p  above  are  the  same  as  in  the  ordinary  blast 
wave  theory,  the  expression  for  is  obtained  so  as  to  fit  the  boundary  condi¬ 
tion  to  B-  in  Eq.  (  3.  32) ,  which  can  be  expressed  as 
6 


Q  _  nL- 
eo  "  2irR  » 


I:  Electric  current 


and  accordingly  we  have, 


b(l)  =1  . 


Substituting  Eqs.  (3.  33)  into  Eqs.  (  3.18) ,  (  3. 19)  and  (  3.  20) ,  we  get  equations 


similar  to  Eqs.  (  2.  25) : 


h(-|f  +  (f-x)^)  = 


5.  k 


V  dx 


X  dx 


(xb) 


j  ,  .  ,  dh  u  #  df  .  f  V 


-Kg  +  (f-x)^+vg(§  +  ^)  =V(V-l)6R^[^^(xb)]^  , 


(3.  34) 


where  we  have  defined 


6  = 


JiJL 


2  .,2„2  »  “m  ■  o■^lRU  * 


(3.35) 


4Tr“  PqR“U 

Eliminating  E  from  Eqs.  (  3.  21)  and  (  3.  23) ,  we  get  an  equation  for  b 
z 

R  -~(xb)]  =’;^[(f-x)b]  . 

m  dx‘-x  dx  •*  dx  '■  ■* 


( 3.  36) 


It  is  readily  seen  that  these  equations  (  3.  34)  and  (  3.  36)  are  consistent  with  the 
boundary  conditions  (  2.  27)  based  on  the  strong  shock  assumption  if  RU  =  const. 
(  =  A  ,  say) ,  from  which  we  get 


X=fi^  =  2,  6, 

y  dR  *  ’  m 
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It  turns  out  that  the  system  of  equations  have  two  intermediate  integrals;  one 
follows  directly  from  Eq.  ( 3.  ^6)  and  the  other  one  from  the  energy  consideration. 
These  are 

r“  --^(xb)  =  (f-x)b  +  K  , 
m  X  dx  ’ 

-x^  [  T^-..  +  xf[-^,  +-^h  -  6R“  b^(bx)  =  K'  , 

yi  Y”l)  2  •*  *■  Y“1  2  ■*  m  dx  * 

where  K,  K'  are  the  integration  constants.  For  further  developments  in  seeking 
the  solution  of  the  system,  it  was  found  that  the  four  conditions  given  in  Eqs. 

(  2.  27)  and  b(  1)  =1  are  not  enough.  An  additional  condition  [xb(x)  =  1 
(a  concentrated  imit  axial  current)  must  be  satisfied  to  have  a  similarity  solution 
which  is  finite  and  which  applies  to  the  entire  domain  0  <  x  <  1  . 

It  is  interesting  to  utilize  the  solution  to  check  the  validity  of  the  approxi¬ 
mate  solution  given  by  the  method  considered  above  for  small  cr  .  It  was  found 
that  both  agree  well  for  small  cr  ,  except  for  the  region  near  x  =  0,  where  the 
purterbation  for  small  o-  was  found  to  be  not  valid. 

3. 4. 2  Magnetohydrodynamic  Blast  Waves 

Electrical  conductivity  is  assumed  infinite  in  this  section.  This  limiting 
case  simplifies  the  problem  to  approximate  the  situation  of  high  conductivity. 

The  interaction  between  the  magnetic  and  the  flow  fields  is  strongest  in  this  case 
and  the  magnetic  field  is  "firozen"  in  the  fluid.  In  consequence  of  this  fact,  the 
fluid  stays  without  magnetic  field  all  the  time  if  no  magnetic  field  permeated 
initially,  and  the  fluid  behaves  as  ordinary  non-conducting  fluid.  The  effect  of 
the  outside  magnetic  field  comes  only  through  the  bounding  surface  acting  as 
magnetic  pressure.  An  example  of  this  type  of  flow  associated  with  the  blast 
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wave  theory  is  seen  in  the  analysis  on  the  transient  state  of  plasma  pinch  given 
by  Kuwabara  (1958,  1963) .  In  this  example  the  plasma  is  assumed  to  be  a  perfect 
conductor  and  there  is  no  magnetic  field  permeated  initially  in  the  plasma  column. 
Two  different  configurations  of  self-and  induced-pinch  are  considered  simultan¬ 
eously.  These  are  caused  either  by  an  azimuthal  or  axial  magnetic  field,  but 
they  act  only  as  magnetic  pressure  to  press  the  plasma  column  inwards  at  ^  =^2 
(  c.  f.  Figure  17) .  This  results  in  a  cylindrical  shock  wave  at  r  =  which 


Figure  17 

Cross  section  of  the  plasma  cylinder 

converges  to  the  center  axis.  The  raotion  of  the  plasma  is  purely  hydrodynamic 
as  described  by  Eqs.  (  2. 1) ,  (2.  2)  and  (  2.  3) .  Assuming  the  strong  shock  condi¬ 
tion  (  2.  27)  and  the  similarity  solution  given  by  Eq.  (  2.  25)  (  in  the  range  of 
X  >  1,  different  from  the  blast  wave  case  of  0  <  x  <  1) ,  Kuwabara  examined  the 
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sheet  current  at  r  1^2  ^^^cted  from  the  similarity  configuration.  He  also 
examined  the  converging  feature  of  the  column,  which  revealed  that  the  amount 

small  during  the  process  of  pinch.  This  fact  is  consistent  with 
the  assumption  considered  by  Allen  (1957)  in  which  he  analized  the  problem  by 
use  of  the  snow  plough  model. 

With  magnetic  field  permeated  initially,  the  fundamental  equations  for  the 
cylindrical  type  flow  given  in  Eqs.  (  3.18)  -(  3.  24)  are  simplified  under  the  as¬ 
sumption  of  infinite  conductivity  and  we  get  from  Eqs.  (  3.  20)  -(  3.  24) , 


)  3 

<  ir  +  3r“»e=°*  . 


1  a 


(3.  37) 


while  Eqs.  (3.18)  and  (3.19)  are  unchanged. 

Two  different  cases  must  be  distinguished  depending  on  the  condition  at  the 
shock  front.  First,  we  consider  the  case  of  o-  =  <»  inside  the  blast  wave  but 
O'  =  0  outside,  which  may  correspond  to  the  case  of  ionizing  shock  front.  Since 
the  magnetic  and  the  electric  field  should  remain  continuous  in  this  case,  we  have 
the  ordinally  shock  condition  ( 2.6)  to  give  the  jumps  in  the  velocity,  the  pressure 
and  the  density.  The  second  case  assumes  a  =  <»  everywhere,  where  upon  the 
magnetic  field  may  be  discontinuous  at  the  shock  front  resulting  from  a  sheet 
current  there.  More  details  of  the  shock  condition  associated  with  this  case  may 
be  found  from  the  following  conservation  relations  of  mass,  momentum  and  energy, 
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< 


(p(U-u)]^^^=  PO  U  , 

r  2  12,  2  2  12 

[p  +  p(U-u)  +-B  =  +-B„  , 

[  p(  u  -  U)  {| ( u  -  U)  ="  + f  )  +  ^  b2(  u  -  U)] 


-Bq  u  , 
JJL  0  ' 


supplemented  by  the  "frozen"  condition, 


[p/B]^^R  =  Pq/^o 


The  conditions  above  lead  to  the  following  conditions 


<  P*>  r=R/’o  "  ^  “Jr=R/‘  ^  ' 

Po/p>r=R  =  3o/(B),,R=l-(u],^^/TJ  , 


(3.  38) 


where  we  have  defined 


jjc  12 

P  =P+^B  ,  y  .C 


"‘/ U^,  C*  =  Y  p/p,  A  =  U-ZT^^Bq  . 

It  is  interesting  to  notice  that  the  condition  (3.  38)  in  the  special  case  of  y  =  2  , 


we  have 


^(u)r,R=fu(l  +  y")  , 


.  1  ❖  *-l 


^  P^r=R  ^  2y‘)”^  , 


(3.39) 


which  has  the  same  form  as  the  ordinary  shock  condition  given  in  Eq.  (  2,  6)  for 


❖ 


Y  =  2,  provided  p  is  replaced  by  p  .It  should  be  recalled  that  y  =  2  is 
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known  to  approximate  the  plasma  state  in  a  uniform  magnetic  field  (  c.  f.  Section 
2.  3.  2) . 

Among  the  various  specialized  configurations  associated  with  applications, 
we  first  consider  the  case, 

B^Q  =  Mnst.,  Bjo=0,  Eq  =  0 

which  corresponds  to  the  situation  above  in  the  problem  of  estimating  the  effects 
of  applied  axial  magnetic  field  on  the  blast  wave  from  exploding  wire  ( Sakurai, 
1962b).  O’  is  assumed  infinite  everywhere  and  the  shock  condition  is  given  by 


Eq.  (  3.  38) .  Equations  (  3.18) ,  (  3. 19)  and  (  3.  20)  are  simplified  in  this  present 
configuration  to  give. 


Du  iJ./ 

Dt  ”  “  p  ar  ‘*2ir  z' 


p  Dt  “  YP  Ert  ~  B 


(3.40) 


The  problem  is  simplified  with  a  further  assumption  of  y  ~  2 
become  « 

I^__L-P2 _ 

Dt~  pat’pDt'  ^Dt"  r8r''  ’ 


and  Eqs.  (3.  40) 
(3.41) 


The  system  of  equations  (  3.  41)  is  the  same  as  that  given  by  Eqs.  (2.1),  (2. 2) 
and  (  2.  3)  for  ordinary  non-conducting  fluid  with  a  =  1,  y  -  Z  ,  Thus  this  sy^stem 
supplemented  by  the  condition  (3.  39)  provides  the  same  solution  as  considered 
in  Chapter  2  on  the  blast  wave.  So  we  can  get  the  cylindrical  blast  wave  solu- 
tion  under  the  axial  magnetic  field  by  simply  replacing  p  by  p  in  the  solution 
given  in  Chapter  2.  It  turns  out  that  Eq.  (  3.  41)  with  Eq.  (  3.  39)  represents  a 
more  general  motion  of  the  c;ylindrical  shock  wave  of  a  progressive  type  (Courant 
and  Friedrichs,  1948)  and  ail  solutions  given  for  the  ordinary  gas  are  valid  in  this 


78- 


case.  In  fact  Kuwabara  ( 1962)  used  this  system  to  analyze  the  transient  state 
of  pinch  phenomenon  of  plasma  column  in  which  the  axial  magnetic  field  is  per¬ 
meated  initially.  Many  other  problems  may  be  treated  by  this  scheme  by  choos¬ 
ing  \  ss  2  to  get  a  rough  idea  of  the  features. 

It  is  also  known  in  this  case  of  y  =  2  that  the  more  general  situation  of 

BzO=  const.,  ^  _ 

can  be  represented  by  the  similarity  solution  for  the  strong  shock  assumption 
(Korobeinikov,  1962).  Also  the  similarity  solution  for  general  values  of  y  is 
possible  if  the  magnetic  field  is  purely  azimuthal  such  as  given  by 

BzO  =  B0O‘*=  Bo  =  0  , 

(Pai,  1958,  Cole  and  Greitinger,  1962,  Korobeinikov,  1962). 

There  is  another  interesting  application  of  blast  wave  theory  to  the  phenomenon 
called  ' ’inverse  pinch"  ( Anderson  eU  ^.1958,  Korobeinikov  and  Ryazanov,  1960, 
Cole  and  Greifinger,  1961,  Liepmann  and  Vlases,  1962).  In  this  situation  we  have 
an  applied  axial  magnetic  field  permeated  in  an  ideal  gas  of  infir  .  conductivity, 
and  an  electric  current  through  the  axis  ( wire  does  not  explode) .  Thus  we  get 
an  expanding  cylindrical  vacuum  region  pushed  by  the  magnetic  pressure  pro¬ 
duced  from  an  axial  electric  current  on  the  cylindrical  surface  which  is  returning 
along  the  axis.  The  magnetic  pressure  acts  as  a  cylindrical  piston  and  results 
an  expanding  cylindrical  shock  wave  ahead.  It  is  interesting  to  note  that  there 
enters  no  characteristic  length  or  time  in  the  initial  or  boundary  conditions  of 
this  problem.  Thus  the  solution  should  have  the  property  of  psuedo-steady.  In 
which  case  it  could  be  expressed  by  a  similarity  solution  of  blast  wave  type, 

where  the  shock  velocity  is  constant  and  the  electric  current  is  proportional  to 
the  time. 
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3,  5  Cavitation 


By  way  of  studing  the  cavitation  problem  in  water,  Hunter  ( 1960)  introduced 
a  similarity  solution  in  the  final  stage  of  a  collapsing  empity  spherical  cavity. 

The  similarity  however  requires  that  the  velocity  of  sound  vanishes  at  the  cavity 
wall.  Although  this  requirement  is  satisfied  with  the  fluids  obeying  a  polytropic 
law  of  the  form  P^Cp”'  with  m  the  index  ( Hunter,  1963),  it  is  usually  not  the 
case  in  water  where  the  equation  of  state  is  conveniently  given  by  the  Tait^equa** 
tion, 

B  ‘po  ’ 

where  B  is  a  slowly  varying  function  of  entropy  and  y  a  value  of  7,  say. 
The  velocity  of  sound  c  is  given  by 

dp  entropy  ^ 

and  it  does  not  vanish  at  the  cavity  wall  since  p  remains  finite  as  the  pressure 

p  goes  to  zero  there.  So  the  similarity  solution  satisfies  the  boundary  condition 

only  in  the  limiting  case  of  p  »  B  and  when  the  radius  of  the  bubble  R  becomes 

small  ( or  dR/dt  becomes  large) .  To  improve  the  solution  for  larger  R,  and 

still  take  into  account  the  effect  of  the  finite  density  at  the  cavity  wall,  Holt 

and  Schwarz  ( 1963)  used  a  perturbation  method.  This  method  involved  expanding 

dR  “1 

the  solution  in  a  series  ^ ®  fashion  similar  to  that  given  in  Section 
2.  3.  It  turns  out  that  the  procedure  is  not  straightft  a  because  of  singular¬ 
ities  involve<^  and  must  be  modified  by  introducing  Lighthill*  s  technique  ( 1949) 
to  avoid  the  singularities. 

The  velocity  and  the  velocity  of  sound  are  given  by  the  following  equations 
of  momentum  and  continuity, 
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r8u  8u  ^^. 

dt  ar  Y-1  9r 


I 


dc^  .  dc^  .  .  ..  2,au,^.  ^ 

- +u-r— +(y-1)c  t'gjr  +  7")  = t 


(3.42) 


at 


dr 


where  u,  Cj  r,  t  have  been  expressed  in  dimensionless  forms.  Equations 


(3.42)  are  supplemented  by  the  boundary  conditions, 
u=R,  c  =  l  at  r  =  R  . 


(  3.43) 


It  was  shown  by  Hunter  (1960)  that  Eqs.  (  3.  42)  have  a  similarity  solution  such 
that  I'/k,  cV^^  sire  functions  of  r/R  (  =x)  only,  and  RR/R^  =  K,  (  3.44) 
where  K  is  a  constant  to  be  determined  from  the  regularity  of  the  solution. 

Now  the  solution  u/R  satisfies  the  boundary  condition  u  =  R  by  setting 
(u/R  )  =  1,  while  (c  /R  )  has  to  be  put  zero  since  R-*oo  as  R“*’0  , 

which  is  however  compatible  with  the  condition  c  =  1  only  for  small  R  .  To 
find  the  solution  valid  for  a  wider  range  of  R,  Holt  and  Schwartz  (1963)  introduced 
a  transformation  of  variables  similar  to  that  given  in  Section  2,  2, 


•  2*2 

u  =  Rj^  f ( X,  y) ,  c  =  R^  g(  X,  y)  , 


(c) 


(  3.45) 


X  = 


’Si’'"’  ’Si 


r~R 


’'h  ’ 


where  R^^  signifies  Hunter' s  R  defined  in  Eq.  (  3,  44) .  Equations  (  3.  43)  are 
then  transformed  into. 


df  t,  af  1  ag 

Kf  +  (  f-x  -  Ky  ~  =  - — ;  ^  , 
ay  ay  Y-1  ax  * 


(3.46) 


2Kg  +  (£-x)|^-Ky|^  =  - 


(Y-1)  +  > 


which  are  very  similar  to  Eqs.  (  2. 15)  and  their  solution  for  small  y  may  be 


found  in  a  similar  manner  as  in  Section  2.3. 
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Expanding  f  and  g  in  series  of  y  , 


f  =  X)  +  ^Nx)  +  y¥  X)  +  . . .  , 


g  =  g^°Nx)  +  yV^Nx)  +  y'^g^  ^^x)  +  ...  , 


(3.47) 


we  get  from  Eqs.  (  3.  46)  the  following  systems  of  equations, 


.(  0) 


(3.  48) 


j ,  0)  2.g(  «)  ]  ^  = ,  v-1)  Kf‘  “>  g'  O)  -  Eg'  «>  (f'  “>  -x)  (  K+  ^  f«»  )  . 


^[(f«'>-.)V°>]4^’=Pi(x)g<«4P,(x)f">  , 


< 


dx 

.(1) 


(  3.  49) 


where  ^^(x),  ...,  |3^(x)  are  functions  given  by  f^^\  g^  .  R  is  considered 
as  a  function  of  y  only  and  is  expanded  as 


R  =  Rjj(  1  +  y  +  . . . )  , 


(3.50) 


where  a^,  ...  are  constants  to  be  determined. 

Utilizing  the  expansion  formulae  given  in  Eqs.  (  3.  50)  and  (  3.  47) ,  the  boundary 
condition  in  Eq.  (  3.  43)  are  expanded  as, 

u=R^[f^®Nl)  +y^{f^®^'(l)a^  +  f^^Nl)}+  ...]  =^[l  +  y^(l-2K)+...]  , 


( l/y^)  [ 1)  (i)®^  +  g' ^^( D)  +  . . .]  =i  , 


2  /  ( 0)  • 


( 1) 


from  which  we  get 


f^°Nl)  =  1,  g^°Nl)  =0  ; 

f^^Nl)  +ajf^°^'(l)  =aj(l-2K),  g^^N  1)  +  '  (1) 


(3.51) 

=  1;  (3.52) 
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The  first  approximation  of  the  problem  is  given  by  Eqs.  (  3,  48) ,  (  3,  51) . 

Unlike  the  blast  wave  theory  in  Section  2.  2,  Eqs,  {  3.  48)  have  two  singularities 

resulting  from  zeros  of  the  denominator  ( f^  -x)  ^  ,  which  are  one  at 

X  =  1  and  the  other  one  is  somewhere  at  x  =  x  >  1  (  say)  •  A  singularity  at 

s 

X  =  x^  appears  also  in  the  blast  wave  equation  (  2.  25) ,  but  we  had  no  difficulty 

since  we  were  only  interested  in  the  region  x  <  1  .  It  can  be  shown  that  there 

exists  a  regular  solution  at  x  &  1  for  any  values  of  K  •  However  a  special 

choice  of  K  makes  the  solution  regular  at  x  =  x  ,  and  we  can  determine  the 

s 

value  of  K  as  well  as  the  value  of  by  using  this  fact.  The  procedure  was 

performed  numerically  by  Hunter  ( 1960)  for  y  =  7  and  K,  x  values  thus 

s 

determined  are  -0,801  and  1,  51  respectively.  Note  that  the  procedure  is  similar 

to  that  for  finding  proper  solution  of  some  other  problems  like  "converging 

shock"  ( Guderley,  1942) ,  "shock  wave  at  the  edge  of  a  gas"  (  Gandel*  man  and 

Frank-Kamenetskii,  1956,  Sakurai,  1960) , 

With  the  value  of  K  as  well  as  the  solution  f^  ,  g^  so  determined 

Eqs,  (3,  49) ,  (  3,  52)  provide  the  second  approximation  to  the  problem.  Note 

that  the  factor  ( f^  -x)  ^  -  g^  appears  also  in  Eqs.  (  3.  49)  and  becomes  zero 

at  X  =  1  and  x  =  x^  .  The  equation  is  singular  there.  Since  the  right  hand  sides  of  the 

equations  (3.49 )  do  not  vanish  at  these  points.  This  indicates  that  the  expansion 

solution  in  y  is  not  valid  there.  However,  it  was  found  by  Holt  and  Schwartz  (1963 )  that 

the  regular  solution  near  x  =  x  could  be  obtained  by  modifying  the  expansion 

s 

procedure  with  the  use  of  Lighthill'  s  technique.  This  involved  the  changing  of 

the  independent  variables  from  x,  y  to  z,y  .  Where, 

2  4 

X  =  Xg  +  z  +  y  qj(  z)  +  y  q^C  z)  +  . . .  , 
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and  the  functions  Qp  q^,  • . .  are  to  be  determined. 

The  expansion  in  y  is  then  assumed  near  z  =  0  as 
f  =  z)  +  z)  +  . . . 
g  =  (  z)  +  g^  (  z)  y^  +  . . .  . 

The  equations  for  f^  ,  g^  are  the  same  as  given  in  Eq.  (  3.  48)  in  terms  of 
X  .  However  the  equations  for  f^^^,  g^  are  different  from  Eq.  (3.49)  and 


are  given  by 

r 


=  +  1^2  ■*■  ^2^1  » 


(  3.53) 


where  6p  6^,  6^  are  all  known  functions  defined  by  K 


We  are  required  to  find  regular  expansion  solution  of  Eq.  (  3.  53)  as 

f^  =  Art  +  A,  z  +  . . . ,  g^^^  =Brt  +  B,z  +  ...  ,  (3.54) 


so  that  the  right  hand  side  of  the  equation  (  3.53)  vanishes  at  z  =  0  .  Sub¬ 
stituting  the  expansion  (  3.  54)  into  Eq.  (  3.  53)  and  assuming  that  q^=constant 
so  that  qj  =  0,  we  get  the  following  condition, 


(P^  =  P.(l),  6^  =  6.(1),  i  =  1,...4) 


from  which  we  have 


P2«3-Vl  '¥i-¥3'  ^1^4 -''2  Pi  ’ 


(3.55) 
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where  the  denominators  are  all  proved  to  be  nonzero  at  z  =  0  .  The  higher- 
order  coefficients  A^,  Bj,  ...  can  be  expressed  as  known,  nonzero  multiples 
of  Aq  .  Hence,  expansions  of  the  form  (  3.  54)  can  be  found  in  which  the 
coefficients  depend  on  one  arbitrary  parameter.  In  this  way,  one-parameter 
families  of  regular  integral  curves,  for  f^^^  and  g^^^ ,  can  be  found  in  the 
neighborhood  of  z  =  0  . 

We  can  also  utilize  Lighthill' s  technique  to  get  a  regular  expansion  solu¬ 
tion  near  x  =  1  and  again  this  solution  involves  a  parameter.  The  features, 
however,  are  slightly  different  from  that  for  near  x  =  x^  described  above. 

Assume  the  expansion  near  x  =  1  , 


< 


X  =  1  +  ;  +  y  q^(C)  +  ... 

f  +y^f^^Nc)  +  ... 

g  =  U  +  y^  g^  +  • . .  , 


which  leads  to  equations  for  f^  ,  g^  similar  with  Eqs.  (  3,  53) .  We  expect 


a  regular  expansion  solution  near  ^  -  I  expressed  as 
f^  '  ^0  ^1^  +  . . .  , 

g^  '  ^0  **’  ^1^  +  . . .  > 


(  3.56) 


and  this  requires 


*^1  ^0  ^2  ^0  ®1  ^1  ®  » 


-g^  (1)  Dq  ^4  ^0  H  • 


(  3.57) 

(  3.58) 


We  have  a  relation  between  and  given  by  Eq.  (  3.  52)  , 

n  -  1  •4-  9  .  1 

2K  +  r  '  ( 1)  "  1 


(3.59) 
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which  reduces  Eq.  (  3.  37)  to 


+  2(K+y)]  Dq- 


2K-l+f 


(0) 


Y-1 


ii). 


2K-i 

Y-1 


=  0 


(3.60) 


It  can  be  shown,  however,  that 

Yf^°^*(l)  -  3+  2(K+  y)  =  0 

and  the  condition  (3 . 60)  is  satisfied  for  any  values  of  Dq  if  is  determined 
by 

Y-1  Y*“l  ® 

Since  Cq  is  determined  by  Eq.  (  3. 5  9)  for  a  given  and  Dj^  is  given 

i  t 

by  Eq.  (  3.  58)  and  so  forth  for  Cp  D^,  . . . ,  we  may  have  a  regular  expansion 
for  any  value  of  Dq  • 

The  parameter  Dq  as  well  as  the  other  parameter  Aq  involved  in  the 
other  solution  (  3. 54)  can  be  determined  by  fitting  the  solution  Eq.  (  3. 56)  with 
Eq.  (  3,  54)  originated  from  x  «  x 

Correction  to  the  cavity  wall  velocity  ap  which  appears  in  Eq.  (  3. 50) 
can  be  obtained  from  Eq.  (  3.  52) ,  with  use  of  the  value  of  determined  above. 


3.6  Thunder 

Although  one  may  easily  imagine  that  the  thunder  is  a  kind  of  blast  wave 
generated  by  lightning,  it  was  realized  only  very  recently  that  the  phenomenon 
was  actually  different  from  the  ordinary  sound  wave  in  some  aspects.  The 
attention  to  the  phenomenon  has  hitherto  been  concentrated  mostly  on  the 
electrical  character  of  the  lightning  stroke.  Also  the  way  of  observing  this 
unpredictable  occurance  of  the  phenomenon  has  made  it  difficult  to  study  the 
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exact  feature  of  the  phenomenon. 

Quite  recently  Kitagawa  and  Kobayashi  ( 1963)  examined  their  recorded  data 
on  the  observed  thunder  and  the  change  in  electric  field  which  were  both  pro¬ 
duced  by  the  same  lightning.  This  examination  revealed  that  the  time  interval 
between  the  lightning  and  the  succeeding  thunder  is  in  fact  shorter  than  the  one 
expected  from  the  propagation  by  the  sound  speed.  This  fact  suggests  that 
thunder  is  a  cylindrical  blast  wave  generated  by  the  lightning  as  a  line  source. 
They  also  estimated  the  range  in  which  the  propagation  velocity  was  greater 
than  that  of  sound  and  concluded  that  it  spreads  over  a  few  hundred  meters. 
Encouraged  by  this  fact,  they  tried  to  estimate  the  characteristic  length 
given  in  Eq.  (  2. 14) .  Lightning  usually  consists  of  sequential  strokes  and  if  we 

denote  the  energy  from  one  stroke  as  W  ,  we  get  from  Eq.  (  2. 14) , 

s 

f 

where  L  is  the  typical  length  of  the  stroke.  To  get  W  ,  they  utilized  the 

s 

formula  given  by  Remillard  (1960)  as 

Wg  =  pLK^(^)^/8ir  , 

where  K  is  a  constant  given  as  10  cm  /amp  .  di/dt  is  the  time  rate  of  the 

9  10 

stroke  current  whose  value  ranges  somewhere  in  the  neighborhood  of  10  -^10 
amp/sec  .  They  obtained  value  of  R^  as  0, 1'^lm  .  These  figures  mean 
that  actually  no  damage  may  be  expected  from  the  direct  effect  of  the  thunder  as 
a  blast  wave  except  in  the  near  vicinity  of  the  lightning.  But  it  should  be  re¬ 
called  that  the  range,  where  the  characteristics  of  the  wave  is  different  from 
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that  of  the  sound  wave,  spreads  much  greater  than  the  length  (  See  Sub¬ 
section  3.1.2).  Although  it  is  not  easy  to  state  a  definite  region  for  his,  it 
is  estimated  by  Thome  (1962)  that  the  distance  where  U/c  becomes  1.1  is 
about  100m,  for  the  typical  lightning  stroke.  This  can  be  compared  with  the 
conclusion  about  the  range  given  above  by  Kitagawa  and  Kbbayashi.  This  wide 
range  of  results  leads  us  to  believe  that  the  phenomenon  may  be  more-  appropriatly 
studied  by  taking  into  account  the  non-linear  effects.  It  is  expected  that  this 
approach  will  give  quite  different  features  to  the  situation  than  the  simple 
acoustic  theory. 

Thome  ( 1962)  studied  the  problem  extensively  based  on  blast  wave  theory 
which  is  regarded  as  the  simplest  possible  model  of  the  phenomenon.  It  is 
noted  however  that  the  range  of  y=c  /U  given  in  Eq.  (  2. 10) ,  which  is  im¬ 
portant  to  this  problem  is  different  from  the  one  to  the  general  blast  wave 
problem.  In  this  problem  we  are  interested  in  a  region  where  y  is  almost  1  , 
to  which  the  series  expansion  solution  is  apparently  inadequate,  and  a  solution 
valid  for  a  wider  range  is  needed.  For  this  purpose,  Thome  utilized  the  approxi¬ 
mation  solution  given  in  Section  2.  4  ( Eqs.  (  2.  58)  -(  2.  68) ) .  Another  difficulty 
is  that  the  pressure  wave  form,  needed  in  making  a  comparison  with  the  observa¬ 
tional  data,  is  to  be  described  as  a  function  of  time  by  an  observer  at  a  fixed 
point.  However,  all  results  given  in  the  blast  wave  theory  are  functions  of 
K,  y  and  are  appropriate  to  give  relation  to  fixed  time  observation,  and  the 
transformation  necessary  to  get  the  relation  is  not  straight  forward.  This  dif- 
riculty  can  be  avoided  since  we  are  interested  only  in  a  region  where  y  is 
very  close  to  1  by  making  use  of  the  fact  that  the  formula  can  be  simplified. 
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The  solution  of  Eq.  (  2.  66)  ne='r  y  =  1  is  found  as, 


^^(Y-1)(«+1)  "  ^^2^  ^>  +  3'*’  ^2  .^  a+1^  +0[{l-y)  ] 

With  use  of  this  expression,  we  get  from  Eq.  (  2.  68) , 


-2^^ - 1 - 

^  ^2^^  (a+l)(a+3) 


(1-y)  +Of(l-y)^3 


or 


y(‘^)  ^  ^  ( 1-y)  +o[(i-y)^] 

Y  -1 


for  a  =  1  , 


which  gives 


R~7KC  , 


(3.61) 


where  we  have  defined  KC 


=  4n/2Ir  . 


Y  0 

The  approximate  formula  (3 , 61)  proves  to  be  good  up  to  about  U/c  -^1. 1  (or 
y  0,  826) ,  when  compared  with  the  numerical  solution  of  Eq.  (  2.  66) .  This 
is  accurate  enough  for  the  present  purpose.  As  noticed  in  Section  2. 4,  Eq. 

(  2. 66)  is  not  precise  enough  to  give  the  exact  behavior  of  J  near  y  =  1  . 
Therefore  Eq.  {  3.61)  is  different  from  the  exact  asymptotic  behavior  near  y  =  1  , 

-3/4 

which  should  beasR  «<l-y.  Nevertheless  the  formula  may  provide  some 
rough  picture  as  to  the  characteristics  of  the  non-linear  nature  of  the  phenomenon 
which  differ  from  those  features  given  by  the  sound  wave  theory.  The  procedure 
may  be  improved  by  using  the  exact  asymptotic  expression,  which  requires  us 
to  determine  unknown  quantities  involved  in  the  expression  in  some  way  by 
fitting  it  with  the  known  solution.  It  should  also  be  recalled  that  we  are  mainly 

interested  in  the  pressure  distribution,  to  which  the  approximate  formula  given 

in  Eq.  (  2,  62)  is  rather  good  in  whole  range  of  y  . 
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Now  k  in  Eq.  {  3. 61)  has  the  dimension  of  the  time,  and  it  turns  out  that 
K  is  very  small  compared  with  the  time  of  arrival  for  most  of  the  observation 
made  at  points  reasonably  far  away  from  the  lightning.  Consider  the  case  of 


R  =  0.1  -  Im  given  above.  By  taking  c  “  340  m/sec,  y  =  1.  4,  k  becomes 


10  ^  -  10  ^  sec.  (Thome  (1962)  used  the  energy  value  of  2.1x10^^  joules 

3  "2 

given  in  Schonland  ( 1950)  and  6. 0x10  m  for  L,  which  gives  2x10  sec 
for  k,  which  is  slightly  greater  because  of  the  larger  value  for  energy) .  In 
any  case,  the  observation  time  for  an  observer  at,  say,  100  meters  from  the 
lightning  will  be  about  0.  3  second  which  is  more  than  ten  times  k  •  Utilizing 
this  fact,  the  transformation  of  the  pressure  function  given  in  Eq.  ( 2. 62)  for 
the  fixed  time  scheme,  can  be  performed  for  t »  k  .  The  result  neglecting 
small  order  terms  in  k/t,  is 


.-2 


10 


’  c  t 


which  gives  the  pressure  distribution  observed  at  a  distance  from  the  lightning 
stroke  r  as  a  function  of  the  time  t  •  It  is  more  convenient  to  rewrite  the 
equation  in  the  following  form. 


pr  =  p^r  + 


(V+IXAI)  ^  c 


where  we  have  used  the  expression  for  k  and  given  in  Eqs.  (  3. 61)  and 
(  2.14)  •  This  relation  can  be  expressed  in  a  single  graph  of  pressure  versus 
time  for  all  observed  distances  by  including  r  as  a  scaling  factor  on  the  axis. 
Such  a  curve  given  by  Thome  ( 1962)  for  E  ^  =  3.  5x10^  joules  is  reproduced  in 
Figure  18,  where  r  is  in  kilometers. 

These  results,  as  well  as  the  frequency  spectrum  of  the  pressure  wave  form 
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Time:  (Scale)  (rinkm)  seconds 


given  above  were  compared  with  observational  data  ( which  are  rather  rare)  and 
also  with  the  theoretical  study  based  on  linear  theory  by  Remillard  ( I960) .  It 
was  found  that  the  results  here  are  significantly  different  from  those  given  on 
a  linear  basis,  but  support  recent  observational  data  by  Arabadzhi  ( 1932) ,  while 
there  are  disagreements  with  earlier  observations  (  Schmidt,  1914) . 

It  is  rather  surprising  that  very  few  studies  have  been  undertaken  in  the 
acoustics  of  thunder  from  the  observational  or  the  theoretical  standpoint.  It 
will  be  very  interesting  to  see  further  developments  in  this  fascinating  field  of 
study,  which  seems  to  be  in  its  infancy. 

Acknowledgement  s 

This  report  was  prepared  when  I  stayed  at  the  Mathematics  Research  Center, 
University  of  Wisconsin,  while  I  was  on  leave  of  absence  from  Tokyo  Electrical 
Engineering  College,  and  it  is  a  pleasure  to  acknowledge  the  support  of 
Professors  R.  E.  Langer  and  J.  B.  Rosser,  Directors  of  the  Center,  and  Dr.  Y. 
Niwa,  President  of  the  College.  I  would  also  like  to  thank  Mr.  D.  Percy, 
Assistant  Director  of  the  Center  for  his  kind  arrangements  of  preparing  the 
manuscripts.  I  am  grateful  to  Mr.  F.  R.  Brown,  Chief  of  Nuclear  W^eapons 
Effects  Division,  Waterways  Experiment  Station  for  his  interests  in  this  article, 
and  to  Mr.  J.  Kirkland  in  the  Division  for  reading  and  correcting  the  manuscripts. 

Finally  I  wish  to  express  my  deepest  gratitudes  to  Professor  M.  Holt,  who 
originally  suggested  to  write  this  article  for  a  chapter  in  “Recent  Advances  in 
Fluid  Mechanics"  and  gave  invaluable  help  with  continuous  encouragements. 


-92- 


REFERENCES 


1.  AUen,  J.  E.  (1957),  Proc.  Phvs.  Soc.  B70,  24. 

2.  Anderson,  O.  A.,  H.  P.  Furth,  J.  M.  Stone,  R.  E.  Wright,  (1958).  Phys. 
Fluids.  1,  489. 

3.  Arabadzhi,  V.  I.  ( 1952) .  Dokl.  A)cad.  Nauk.  SSSR.  12,  377. 

4.  Bennett,  F.  D.  (1958a).  Phvs.  Fluids.  1.  347. 

5.  Bennett,  F.  D.  (1958b).  Phvs.  Fluids.  1,  515. 

6.  Bennett,  F.  D.  (1959e).  Phvs.  Fluids.  1,  470. 

7.  Bennett,  F.  D.  (1959b).  In  "Exploding  Wires"  (Chace  and  Moore,  eds.), 

vol.  1,  pp.  211-226. 

8.  Bennett,  F.  D. ,  H.  S.  Burden,  and  D.  D.  Shear  (1962).  Phvs.  Fluids. 

5,  102. 

9.  Bennett,  F.  D.  ( 1962) .  Phvs.  Fluids.  1,  891. 

10.  Bennett,  F.  D.  (1964).  Phys.  Fluids,  to  appeared, 

11.  Berezin,  O.  A.  (1960),  Dokl.  Akad.  Nauk.  SSSR.  133,  2396, 

12.  Brode,  H.  L.  (1959).  Phvs.  Fluids.  2,  2l7, 

13.  Chandrasekhar,  S.  (1939).  "Stellar  Structure",  Chapter  IV.  Dover. 

14.  Chace,  W,  G,  and  H.  K.  Moore  (1959),  "Exploding  Wires",  vol.  1;  (1962). 

Vol.  2|  ( 1964) .  Vol,  3.  Plenum  Press,  New  York. 

15.  Courant,  R.  and  K.  O.  Friedrichs  ( 1948) .  "Supersonic  Flow  and  Shock  Waves". 
Interscience,  New  York, 

16.  Ericsson,  U.  and  K.  Edin  (1960).  Phvs.  Fluids.  893. 

17.  Filter,  W.  S.  (1960).  Phvs.  Fluids.  _3>  444. 

18.  Friedman,  M.  P.  (1961),  J.  Fluid  Mech. ,  11,  1. 


-93- 


19.  Gandel*man,  G.  M.  and  D.  A.  Frank-Kamentskii  (1956).  Dokl.  Akad. 
Nauk.  SSSR.  107.  8U. 

20.  Gerber,  N.  and  J.  M.  Bartos  ( 1961) .  BRL  Memorandum  Report  1376. 

21.  Glasstone,  S. ,  Editor  (1962).  '*The  Effects  of  Nuclear  Weapons."  United 
States  Atomic  Energy  Commission. 

22.  Goldstine,  H.  and  J.  von  Neumann  (1955).  Comm.  Pure  and  AppI.  Math. 

8,  327. 

23.  Greenspan,  H.  P.  (1962).  Phvs.  Fluids.  ^  255. 

24.  Greifinger,  C.  and  J.  D.  Cole  (1961).  Phvs.  Fluids.  4,  527. 

25.  Greifinger,  C.  and  J.  D.  Cole  (1962).  Phvs.  Fluids.  ^  1597. 

26.  Guderley,  G.  ( 1942) .  Lufo.  19  ^  302, 

27.  Hayes,  W.  (1947).  Quart.  AppI.  Math. .  ^  105. 

28.  Hayes,  W.  and  R.  Probstein  ( 1959) .  "Hypersonic  Flow  Theory.  "  Adacemic 
Press,  New  York. 

29.  Holt,  M.  ( 1956)  •  Proc.  Roy.  Soc.  London,  A234.  89,  111. 

30.  Holt,  M.  and  N.  J.  Schwartz  (1963).  Phvs.  Fluids.  6^  521. 

31.  Hunter,  C.  ( 1960) .  T.  Fluid  Mech.  _8,  241. 

32.  Hunter,  C.  (1963).  1,  Fluid  Mech.  15^  289. 

33.  Jeanmaire,  P.  (1963).  Rep.  No,  9/63.  Emst-Mach-Institut, 

34.  Jones,  C,  W,  ( 1955) .  Proc.  Roy.  Soc.  London.  A228.  82. 

35.  Jones,  D.  L.  (1962).  MBS  Technical  Note  155. 

36.  Jones,  D.  L.  and  R.  M.  Gallet  ( 1962) .  In  "Exploding  Wires,  "  ( Chace  and 

Moore  eds.),  Vol.  2,  pp,  127-144.  Plenum  Press,  New  York. 

37.  Karlikov,  V.  P.  (1955).  Dokl.  Akad.  Nauk  SSSR.  101,  No.  6. 


-94- 


38.  Keller,  J.  B.  ( 1956) .  Quart.  Appl.  Math.  14,  171. 

39.  Kitagawa,  S.  (1963).  JProc.  Int.  Coni.  Atmospheric  and  Space  Electricity. 
Montreux. 

40.  Xochina,  N.  N.  and  N.  S.  Mel'nikova  ( 1958) .  Dokl,  Akad.  Nauk  SSSR. 
122.  192. 

41.  Kochina,  N.  N.  and  N.  S.  Mel*  nikova  ( 1960) .  Akad  Nauk  SSSR.  Prikl. 
Mat,  i  Mek.  24.  213. 

42.  Korobeinikov,  V.  P.  (1958).  Dokl.  Akad.  Nauk.  SSSR  121.  613. 

43.  Korobeinikov,  V.  P.  and  E.  V.  Ryazanov  ( 1959) .  Akad.  Nauk  SSSR.  PrDcl. 
Mat.  i  Mek.  749. 

44.  Korobeinikov,  V.  P.  and  V.  P.  Karlikov  ( i960) ,  Dokl.  Akad.  Nauk  SSSR 
133.  746. 

45.  Korobeinikov,  V.  P.  and  E.  V.  Ryazanov  ( 1960) .  Akad.  Nauk  SSSR.  Frikl. 
Mat,  i  Mek.  24.  111. 

46.  Korobeinikov,  V.  P. ,  N.  S.  Mel'  nikova  and  E.  V.  Ryazanov  ( 1961) . 

"Theory  of  Point  Explosions.  "  ( in  Russian) ,  Gos.  Izd.  Fiz.  -Mat.  Moscow. 

47.  Korobeinikov,  V.  P.  and  E.  V.  Ryazanov  (1962).  Akad.  Nauk  SSSR.  Prikl. 
Mek.  i_  Tec.  Fiz,  4,  47. 

48.  Korobeinikov,  V.  P.  (1962),  ^"Problems  of  Magneto-hydrodynamics  and 
Plasma  dynamics,  "  Riga.,  pp.  197-206. 

49.  Korobeinikov,  V.  P.  and  P.  I.  Chushkui  (1963).  Zhumal  prikladnoi 
Mekahaniki  i  Tekhnicheskoi  fiziki.  4.  48. 

50.  Kuwabara,  S.  i  1958) .  Preprint,  13th  Annual  Meeting,  Phys.  Soc.  Japan, 

III,  2. 


95- 


51.  Kuwabara,  S.  ( 1962) .  Preprint,  17th  Annual  Meeting,  Phys.  Soc.  Japeui, 
V,  30. 

52.  Kuwabara,  S.  ( 1963) .  2?  Phvs.  Soc.  Tapan.  18,  713. 

53.  Landau,  L.  D.  (1957).  In  "Similarity  and  Dimensional  Methods  in 
Mechanics,  "  ( L.  I.  Sedov,  English  translation  by  Hott,  M.  and  M. 
Friedman) ,  Academic  Rress,  New  York,  pp.  296-304. 

54.  Lees,  L.  and  T.  Kubota  ( 1957)  •  J.  Aero.  Sci.  24.  195. 

55.  Liepmann,  H.  W.  and  G.  Vlases  (1961).  Phvs.  Fluids.  4,  927. 

56.  LighthiU,  M.  J.  ( 1949) .  PhU.  Maa.  40,  1179. 

57.  Lin,  S.  C.  ( 1954) .  J.  AppI.  Phvs.  25.  54. 

58.  Morawetz,  C.  S.  (1954).  NYU  Report. 

59.  Ness,  N. ,  J.  B.  Fanucci  and  L.  J.  Kijewski  ( 1963) .  Phvs.  Fluids.  ^ 

1241. 

60.  Oshima,  K.  ( 1960) .  Rep.  Aero.  Res.  Inst.  Univ.  of  Tokyo,  No.  358. 

61.  Oshima,  K.  (1962).  In  "Exploding  Wires"  (Chace  and  Moore,  eds.) 

Vol.  2,  pp.  159-174.  Plenum  Press,  New  York. 

62.  Pal,  S.  I.  (1958).  Air  Force  Office  Scientific  Research,  TN58-215. 

63.  Remillard,  W.  J.  ( I960) .  Office  of  Naval  Research,  Tech.  Memorandum 
No.  44,  NR*384-903. 

64.  Rogers,  M.  H.  ( 1956) .  Proc.  Rov.  Soc.  London.  A235.  120. 

65.  Rouse,  C.  A.  (1959).  ^n  "Exploding  Wires",  Vol.  1.  (Chace  and  Moore, 
eds),  pp.  227-263.  Plenum  Press,  New  York. 

66.  Sakurai,  A.  ( 1953) .  J.  Phvs.  Soc.  Tapan,  8^  662. 

67.  Sakurai,  A.  ( 1954) .  J.  Phvs.  Soc.  Tapan.  ^  256. 


-96- 


68.  Sakurai,  A.  (1955a).  X?  Phvs.  Soc.  Tapan.  10.  827. 

69.  Sakurai,  A.  ( 1955b) .  J.  Phvs.  Soc.  Tapan.  10.  1018. 

70.  Sakurai,  A.  ( 1956) .  i.  Fluid  Mech.  1,  436. 

71.  Sakurai,  A.  (1959).  In  "Exploding  Wires",  Vol.  1.  (Chace  and  Moore, 
eds.),  pp.  264-270.  Plenum  Press,  New  York. 

72.  Sakurai,  A.  (1960).  Comm.  Pure  and  Appl.  Math. .  13.  353. 

73.  Sakurai,  A.  (1962a).  Proc.  12th  Japan  Nat.  Congress  Appl.  Mech.,  157. 

74.  Sakurai,  A.  ( 1962b) .  J.  Phvs.  Soc.  Tapan.  17.  1663. 

75.  Sakurai,  A.  and  T.  Takao(1963).  Rep.  Aero.  Res.  Inst,  Univ.  of  Tokyo, 
No.  3,  A-80. 

76.  Sakurai,  A.  and  T.  Takao(1964).  In  "Exploding  Wires*'.  Vol.  3.  Chace 
and  Moore,  eds. ) 

77.  Schmidt,  (1914).  Met.  Zeits,  487. 

78.  Schonland,  B.  F.  J.  (1950).  In  "The  Flight  of  Thunderbolts",  p.  99. 

Oxford  Univ.  Press,  London. 

79.  Sedov,  L.  I.  (1946).  C.  R.  Ac?d.  Sci. .  U.R.S.  S.  1,  17. 

80.  Sedov,  L.  I.  (1957).  "Similarity  and  Dimensional  Methods  in  Mechanics". 
Translation  in  English  by  Holt,  M.  and  M.  Friedman  ( 1959) .  Academic 
Press,  New  York. 

81.  Spitzer,  L.  (1956).  "Physics  of  Fully  Ionized  Gases",  Interscience, 

New  York. 

82.  Swigart,  R.  J.  ( I960) .  Lockheed  Aircraft  Corporation,  Technical  Report, 
LMSD- 288210. 

83.  Taylor,  G.  I.  (1950).  Proc.  Rov.  Soc. .  London.  A201.  159. 


-97- 


84.  Thome,  G.  D.  (1962).  Thesis  Cornell  Uni v. 

85.  Whitham,  G.  B.  (1950),  Prop.  Rov.  Soc.  London,  A203.  571. 


-98- 


};;•  cr  w 

c  0)  3 

CO  CO  M. 


§■  % 


6!  z 

(D  (D 

3 

5-  3 

(D  a 

i  0) 

o 

a 


o 

00  Q 

•o  ^ 

•a  w 

•  o 


^  §■ 

t-i- 

QJ 

‘<  >3 

0  3 

C  X 
3. 


M-  O 
CO  p 

Cl  p 

<  ^ 
I  3^ 
O  R 
*0 


9  s 

O  50 

s  S 

3  JO 

? 

•  2: 

•*  o 

2  * 

>  ^ 

»L  'O 

i 

O 

i>i  ft 
NJ 


o  f' 

2  3  X 

«  o  S 
o  3  3 

•-^  3  CO 

3^  o  o 

got? 

0^0 

3-  ?  3 

nf 

rf  ^ 

o  3^ 

;? 

0) 


12  >  y 

5.  *0  a 

a  {± 

2  X  Si. 

>  o  o 

Q  o  OJ  m: 

9  jr  3  t!? 

^  3  a  r*’ 

Jl  S-  ^  w 
o  o  £  .Q 

'*  a  OJ 

S  W 

CO  « 


tr  CO 
O  »—  M- 

c  &  3 

CO  CO  13. 


“o  -- 

o  ^ 

§■? 

s-s 

3  *0 

CO  3 

e{  £ 

n>  n> 

3 


fi{  S' 


so  ^ 

*0 
►0  3 

t3  CO 


*-.•  o 

CO  3 

^  3 

o  R 

-D  O' 

n  O 

o.  t-h 


30^ 

m  P  O 

R  f=i-‘  p* 


o  3: 

(D 


§  n 
8  I 

51 

O 

*  2! 
*A  o 

> 

I  s£) 

^  ~4 
H-  ^ 

o  > 

I 

o 

?o 

O 

I 

N 

o 

u» 

s3 


(I)  • 

8  §•  » 

°  I  a 

to  3 

”  2  W 
H-  O  X 

P  “  *0 

M  -  2 


DO  2 
hh  q, 

S' 

t!  3 

^  o> 

> 

S  S 
3  « 

V-M 

«  O 

O  o» 

s  s 


Hj  > 

jr  t3 

E  *0 


d  § 

»-•  55 


w 


O  ji. 


*0  ^ 

II 

CD  <^ 


SI  2: 

CO  o 

•-«  3 
<  O' 
S 

CO  a 

^  p 
p.  -O' 

•  t— ' 
»-•- 

o 

9i 

s£)  Jj. 
00  5 
>0  3 
•o  w 

•  o 


to 

3 

M* 

0 

§ 

Q 

oT  > 

Ql 

JO 

►0  3 

•?•§ 

W  '2 

0 

rt 

0 

?? 

1 

S 

0  3 

• 

•  • 

z 

c  S 

d 

0 

• 

n  ^ 

> 

1 

0  0) 

3  rf- 

vO 

-si 

'i'  w*. 

S  r+ 

p 

o  § 

*a  ^ 


3  o'  CT 
o  9 


a  y 


t-t  ^ 

JJ'  0 


o  f- 

^  (0 

£{ 


g.  Pi 

S  a. 

*>  A 

_ 

>  (0  •  O 

O  O  OJ  Jll 
o  P  3 
C  &  D-  • 

r+  M-  ij-*  M 
O  A 

^  W  *9  C 
to  %#  3  Qj 

s  ^ 

CO  • 


H.  O' 

o  ST 

to  !-♦■ 
•o  g 

o  5 

p  5 

o 


•<  *3 

to  "O 

a  9 

C  A 


o 

Bi 

00  Q 

►o  3 
X5  W 

*  O 


o  3^ 
o  R 

*0  O- 


o  tr 

3  O 


O 

<  ‘t’ 

SI 


a  2 
o  3 
o  n 

0)  PJ 
O  (D 
rf  « 

2?  2 
o  3- 
••  ss: 

2  ° 

i- 

r 

rs»  > 
ts)  d 

I 

o 

JO 

d 

I 

cv 

o 

vO 


DO  >> 

5  I 

£3  I 

<  3 

$  2 
^  n 

CO 

•-3  a 

a  o 

trj  CO 

O  O 
w  D) 

P 


&-  '3  5. 
D.  3  r- 
,  O  Q) 
X  i— 

>•  CD  .  d 

O  R  §  & 

s§  l« 

n' 

p  S  .g  C 

to  V.  93 
3 
3 

w  o 

_ il-3_. 


